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ABSTRACT 


Vibration control is a very important issue in satellites. The new high-resolution digital 
imaging devices are especially sensitive to vibrations. Antennas used in laser communications also 
require a very quiet environment so that their performance is not degraded. The Stewart platform 
is capable of isolating an optical payload from the noisy spacecraft bus. Until recently, only pas¬ 
sive methods were used in all vibration isolation applications. Recent advances in Digital Signal 
Processing techniques made the development of vibration control algorithms possible, but these 
usually require large computational power. This work explores using a computationally efficient 
vibration-isolation method for optical payloads by using hexapods. The method suppresses the vi¬ 
bration at the assigned frequencies and does not affect unassigned frequencies if the plant is linear. 
The mathematical analysis includes convergence analysis and the effect of unassigned frequencies 
in the output. The computational requirements of the algorithm is evaluated and is compared to the 
Multiple-Error Least Mean Square. The method is very robust to nonlinearities; its performance is 
comparable to the Multiple-Error Least Mean Square with a fraction of the computational time and 
memory requirements. Also it requires very little plant knowledge. Theoretical results are verified 
through simulations using a Single-Input/Single-Output plant and a nonlinear hexapod model. The 
controller was also experimentally validated in two different hexapods and the performance was 
found to be similar to or better than the performance obtained with the Multiple-Error Least Mean 
Square method when a noisy reference signal is used. 
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I. INTRODUCTION 


A. THE VIBRATION ISOLATION PROBLEM 

Although the vibration control problem has been studied for a long time, it received increas¬ 
ing attention after the advent of submarines. The main purpose of a submarine is to travel unnoticed 
and the noise generated by its several subsystems can defeat its purpose, making its localization and 
identification easier. As a result, the study of noise and vibration received much attention, especially 
during the Cold War when several techniques for noise and vibration isolation were developed. 

During the sixties, space started being considered as a very strategic area and imaging 
was among the many applications that found their way to the military spacecraft. During the first 
decades of space imaging, vibration was not as important as it is today because of the limitations 
of the imaging sensors available at that time. Nowadays satellites are larger and have much more 
rotating machinery than in recent past. Among the vibration-generating devices present in the recent 
satellites are large appendages, which snap when entering or emerging from an eclipse, solar panels 
that rotate to track the sun, cryocoolers, pumps, reaction wheels, control momentum gyros (CMG), 
reaction jets and magnetic torquers. The vibrations produced by these elements, if transmitted to 
the imaging device, can degrade its performance appreciably. Figure 1.1 shows an example of such 
degradation. 

At the same time, the electronic industry has made astonishing advances in imaging de¬ 
vices. Resolution is now much higher than it was a few years ago, and new optical techniques 
produce lenses and mirrors with less distortion, generating sharper images. Even when distortions 
are present, the images can be post-processed on the ground, using the latest advances in image- 
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(a) Jitter free (b) Jitter present (l/jrad) 

Figure 1.1. Jitter Effeet on an Imaging Deviee (From [1]). 


proeessing teehniques to enhanee them. One famous example is the spherieal aberration on the 
Hubble teleseope’s original primary mirror. Using image proeessing teehniques, the teleseope eould 
still be used while the mirror was being upgraded[6]. 

The most common cause for vibration-induced image degradation is due to vibrations pro¬ 
duced elsewhere and transmitted through the satellite bus to the imaging device. Isolating the noisy 
environment, dampening the structure or isolating the instrument itself by using an isolating mount 
(see Figure 1.2) can reduce this type of vibration. The vibration can also be produced on the quiet 
side when moving parts, such as cryocoolers, must be mounted as closely as possible to the sensor. 

Passive methods are used in most vibration-control and vibration-isolation implementations 
for several reasons: they are robust; they do not consume electrical power; they are very reliable 
and the technology is mature. Passive approaches can be applied to mounts that prevent the vibra- 
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Figure 1.2. The Vibration Isolation Problem (From [1]). 

tions to propagate to or from the spaeeeraft bus or to inerease the dampening of a larger sttueture. 
Unfortunately, passive approaehes do not solve all vibration problems. 

When a vibration isolation mount is used, the vibration has two main paths to migrate from 
the noisy to the quiet side: the vibration isolation medium and the umbilieal. Even though it is 
possible to design very soft passive mounts to preventing the vibration to propagate through the 
mount, prevent the vibration from propagating through the umbilieal is impossible. In the same 
way, passive isolation approaehes are ineapable of suppressing the vibration generated on the quiet 
side. This is espeeially important for infrared imaging for whieh the sensor must be eooled and a 
eryoeooler must be plaeed as elosely as possible to the sensor in order to inerease its performanee. 
Passive solutions ean only perform vibration eaneellation in eertain eases. 

Disturbanee isolation, also ealled “vibration suppression,” is the suppression of periodie 
disturbanees. These disturbanees are usually generated by rotating maehinery and appear as lines 
in the frequeney speetrum. These lines have a power eontent mueh higher than the baekground 
vibration noise. If the frequeney is known and fixed, TMDs (tuned-mass deviees) ean be used, but 
their plaeement on the strueture is eritieal. Thus TMDs are not usually suited for isolation mounts. 
Aetive methods ean deal with frequeneies not known a priory or with frequeneies fluetuating, and 
the plaeement of the aetuators is not as eritieal as when using passive methods. 
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B. ACTIVE VIBRATION SUPPRESSION 


Since meeting the requirements using passive-only solutions is becoming increasingly dif¬ 
ficult, active solutions have been pursued. Current research indicates that the most viable solution is 
the use of a generic local vibration isolator that combines both passive and active isolation. Among 
the several approaches tested, the most promising is the Stewart Platform [7], which can gener¬ 
ate force and torque in any direction. It also exhibits good stiffness properties [2] with actuators 
much smaller than those needed with serial manipulators. The need for localized vibration isolation 
prompted the development of the Vibration Isolation and Suppression System (VISS) [8, 9, 3] and 
the Satellite Ultraquiet Isolation Technology Experiment (SUITE) [10]. 

A hexapod has complex dynamics and a detailed model can be quite complicated. Even 
if a perfect model is developed, space missions usually do not allow the opportunity for scheduled 
maintenance in order to repair aging or failed systems and thus system identification may be needed. 
Eor this reason, all vibration control work found in the reviewed literature uses some sort of adaptive 
method. This approach makes it possible to compensate for model imperfections and aging. 

Since the hexapods are multiple-input/multiple-output (MIMO) systems strongly coupled 
[11], expecting that the best performance can be obtained by using MIMO controllers is reasonable. 
Some methods are available for MIMO plants, such as Multiple-Error EMS[12], Clear Box [13, 4] 
and Repetitive Control [14, 15]. Unfortunately, the computational resources available in satellites, 
limit the use of these methods for vibration isolation of optical payloads. 

One important characteristic common to all hexapods used for vibration isolation as de¬ 
scribed in the literature is the use of smart struts. All of them have six identical struts with an 
actuator and a vibration sensor embedded (geophone or accelerometer), meaning that using one 
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independent single-input/single-output (SISO) eontroller for eaeh strut is possible. This approaeh, 


while simple, has limitations for a system as strongly eoupled as a hexapod, and mixed results have 
been aehieved [11, 16]. 

Eliminating the requirement of having a MIMO eontroller does allow more ehoiees sinee 
several vibration-suppression algorithms were never extended to the MIMO ease. In 2000, Kuo and 
Morgan [17] presented an interesting review of some of the vibration isolation algorithms available 
and, then elassilied them into “referenee-based” and “synthetie referenee.” The first group uses a 
referenee signal eorrelated with the disturbanee while the seeond one uses some other information 
to generate an internal referenee signal. The “synthetie referenee” group eontains is one set of 
eontrollers ealled “noteh filters.” 

Noteh filters have several properties that are very attraetive for vibration isolation in satel¬ 
lites. First, no external signal is used as a referenee. Sinee extra information about the frequeney 
is given or generated on-the-fly, the frequeney information ean be used to generate a synthetie ref¬ 
erenee. This eliminates the requirement of having an additional sensor, at the expense of requiring 
additional information to generate the referenee signal. Another eharaeteristie of this elass of eon¬ 
trollers is that exeellent results ean be aehieved for suppressing vibrations generated by rotating 
maehinery, also ealled tonal disturbanees. 

Among the eontrollers found in the reviewed literature, the one used in 1998 by Bertran 
and Montoro in [18] (based on the work presented by Kuo and Min in 1995 [19]) required the least 
amount of eomputational effort is 

This eontroller is quite simple and does not use the model of the plant explieitly. Although 
this eontroller is very simple, the reviewed literature did not mention this eontroller being applied to 
MIMO systems, espeeially on hexapods. In addition, no extension of Montoro and Bertran’s work 
for MIMO plants was found. More important, no stability analysis was found in the literature that 
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considers a generic plant. However, such analysis is essential to meet the high-reliahility require¬ 


ments of space applications. Another drawback is that the method requires the knowledge of the 
frequency to be suppressed. Other frequencies could be present on the error signal and some inter¬ 
ference can occur. This possibility must be addressed if this controller is to be implemented in an 
actual space system. However, this problem has not yet been addressed in the literature reviewed. 

Although much effort has been spent on developing vibration control and isolation methods, 
our examination of the literature revealed that no one has evaluated the computational effort that 
each method requires. Consequently, a cost analysis must be developed and the vibration-isolation 
method selected must be compared to the Multiple-Error LMS implementation. 

Finally, most modeling work done on hexapods was from a robotics perspective, using 
direct and inverse kinematics. Li and Salcudean presented a model for a hexapod, but assumed that 
the dynamics of the actuators are negligible [20]. Lebret, Liu and Lewis derived a very complete 
model for a hexapod starting from the energy approach, but unfortunately the geometry assumed in 
their work does not apply to the hexapods available at the Satellite Research and Design Center. 
each pair of actuators was connected to the top plate through a single joint. In 1991, Fujimoto, 
Kinoshita et al. derived a model from the robotics perspective [21], but assumed that “the joint does 
not move in high speed,” preventing its application from being used for vibration-isolation purposes. 
A complete model of the hexapod suitable for simulations is needed in order to study the effects of 
the multiple-input/multiple-output (MIMO) nature of the hexapod without the interference of actual 
hardware imperfections. 
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C. PROBLEM STATEMENT 


The goal of this research is to study a computationally efficient algorithm that meets the 
vibration-isolation requirements of optical payloads using hexapods. The method needs to be tol¬ 
erant to nonlinearities, must require significantly less computational resources than the Multiple- 
Error LMS implementation, and the method must perform comparably to the Multiple-Error EMS 
method. Additionally, the method must be able to deal with multiple frequencies. 

A stability analysis and a study of the effect of uncontrolled frequencies must be performed. 
The results must be validated on simulations and experimentally. 

Einally, the computational cost of both the proposed method and the Multiple-Error EMS 
approach must be derived and must be compared in order to verify if the proposed method 
achieves the desired goal of significantly reducing the computational cost over the reference method 
(Multiple-Error EMS). 

D. OVERVIEW 

Chapter I presents the vibration-isolation problem, reviews the literature and states the re¬ 
search objectives. 

Chapter II describes the two different hexapods used to perform the experiments. 

Chapter III provides a detailed state-space model of the Precision Pointing Hexapod. This 
model allows the designer to derive some important properties of the plant and is used to implement 
a detailed simulation program. 

Chapter IV presents the Adaptive Disturbance Canceller method, derives the stability cri¬ 
terion, and studies the method’s performance when multiple frequencies are present. Simulations 
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are used to validate the mathematieal results using both a SISO plant and the model developed in 


Chapter III. 

Chapter V derives and eompares the eomputational requirements for both the Adaptive 
Disturbanee Caneeller and the Multiple-Error LMS approaeh, as well as derives and eompares the 
eorresponding memory requirements. 

Chapter VI validates the results experimentally by using the two hexapods available at the 
Satellite Research and Design Center, at the Naval Postgraduate Sehool. A eomparison between 
the Adaptive Disturbanee Caneeller and Multiple-Error EMS is also presented. 

Einally, Chapter VII presents a summary of the work and eonelusions and some suggestions 


to overeome the main weaknesses of the Adaptive Disturbanee Caneeller. 



II. HARDWARE DESCRIPTION 


A. OVERVIEW 

The Stewart platform, also known as a “parallel manipulator,” was introdueed by Stewart 
in 1965[7] and is eomposed of two plates eonneeted by some (passive or aetive) links. One of the 
plates is eonsidered to be the referenee (base) and the other (top) is the one for whieh the attitude 
and position are to be eontrolled. A diagram of sueh a platform is shown in Figure 2.1. The platform 
ean be eontrolled by ehanging the length of the links. It ean be shown that in order to get six degrees 
of freedom, one needs at least six aetuated struts. One ean also verify that more than six aetuators 
will not inerease the number of the DOF of the manipulator (differently from the serial-linked ease). 
The hexapods have several interesting eharaeteristies, namely, high stiffness and redueed size and 
weight when eompared to serial-link manipulators. A eomprehensive study of hexapods and their 
eharaeteristies ean be found in [2, 22, 23, 24, 25, 26]. 



Figure 2.1. Hexapod View in 3D (From[2]). 
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Section B describes the first hexapod available at the Satellite Research and Design Center, 


the Ultra-Quiet Platform. This hexapod is equipped with high bandwidth piezoelectric and sensitive 
geophone sensors on each strut. 

Section C details the Precision Pointing Hexapod, which has long stroke actuators and is 
capable of pointing and can simultaneously isolate vibration. Since this hexapod is modeled in 
Chapter III, it is described in more detail. 

B. ULTRA-QUIET PLATFORM 

This Ultra-Quiet Platform, shown in Figure 2.2, was conceived for vibration-isolation pur¬ 
poses. It follows a cubic configuration (see Figure 2.3) in order to minimize the cross-coupling 
among the struts [27]. The hexapod rests on a satellite bus mockup (Figure 2.4) where the distur¬ 
bance source is mounted. The whole structure is mounted, with rubber vibration suppressors, on a 
Newport RS4000 floating table in order to provide isolation from the ambient vibration noise. 



Figure 2.2. Ultra-Quiet Platform. 
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Figure 2.3. Cubic Configuration (From [3]). 

The Ultra-Quiet Platform struts are equipped with piezoceramic stack actuators (PZT). 
These actuators have a very large bandwidth (resonance peak around l.lKHz) and a displace¬ 
ment of ±50qm, which is adequate for vibration control. The power to the actuators is provided 
by a PCB Piezotronics 790A06 amplifier (±200F at ± 100mA). Each smart strut is also equipped 
with a Geospace GS-llD geophone, which measures linear velocity. These sensors have a natural 
frequency of l4Hz with a damping ratio of 0.8. The signal conditioning is provided by the CSA 
Engineering Active Vibration Control System (AVCS). The fixtures have passive isolation built-in, 
which provides damping in series with the actuators. 

The disturbance source produces acceleration parallel to the z-axis, as defined in Figure 2.1. 
It uses an Aura Bass Shaker (AST-lB-4, 251T), and it is mounted on top of the mockup, below the 
bottom plate of the hexapod. The shaker signal is produced by a DSP board and is amplified by a 
DC power amplifier {Kepco BOP 20-10M). 

The implementation of the algorithm is made on a dSpace Alpha Combo. This system 
comprises a DS1103 and an Alpha Board. The Alpha Board is powered by a DEC Alpha SQOMHz 
processor with 2MB of RAM. The DS1003 board acts as a slave I/O board and contains a Texas 

11 





Instruments TM320C40 processor running at 50MHz and 512^^; of memory. The setup is shown in 


Figure 2.4. The controller runs on the Alpha Board, while all the I/O functions are executed on the 
DS1003 board. 

All the I/O are performed at 10^ samples/s while the controller and disturbance generator 
run at 10^ samples/s. Since there is not appreciable power contents on the measurements above 
2KHz, anti-aliasing filters were implemented by oversampling the measurements (lOKHz), passing 
them through a Chebychev (type I, corner at 200Hz, 3 order) and finally down-sampling fhe oufpuf 
of fhe fillers lo 10^ samples/s. 

The oulpuls of fhe conlroller are up-sampled lo \0‘^ samples/s and Ihen low-pass fillered 
by a type I Chebychev filler in order lo smoolh Ihe oulpul senl lo Ihe aclualors. The dislurbance 
signal received Ihe same Irealmenl. Filtering, A/D sampling, D/A conversions and communication 
belween Ihe conlroller and Ihe hosl PC occur in Ihe DS1003 board. 

The code is implemented, compiled and downloaded lo bolh boards Ihrough Ihe 
dSpace/Mallab Real-Time Workshop environment The aclual conlroller, as well as Ihe dislurbance 
generator, is custom-coded in C. These lasks are accomplished by using a hosl PC, an Intel-based 
P-III Dell Dimension XPS 600MHz- The same hosl PC is also used to coiled all Ihe relevanl dala 
produced during Ihe experiment The dala are Ihen converted into Mallab’s proprielary formal for 
storage. 
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Figure 2.4. UQP Experimental Setup (From [4]). 
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Figure 2.5. Ultra-Quiet Platform Controller. 
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C. PRECISION POINTING HEXAPOD 


The main limitation of the Ultra-Quiet Platform is its small stroke, whieh prevents it from 
being used for steering. In order to evaluate the possibility of performing both pointing and vibration 
isolation, a hexapod with a longer stroke is needed. To address this need, the Satellite Research and 
Design Center aequired a new hexapod: the Preeision Pointing Hexapod, from CSA Engineering, 
Ine.. The platform is designed to allow both vibration eontrol and pointing. In order to aehieve 
the large stroke required, the Preeision Pointing Hexapod uses voiee-eoil aetuators, whieh provide 
±5mm travel. Aeeelerometers are mounted inline to provide vibration information. The platform, 
as delivered by CSA Engineering, Ine., is depieted in Eigure 2.6. 

This hexapod was seleeted for modeling beeause it has a larger stroke. This hexapod will be 
used for pointing researeh as for well as vibration-isolation researeh and thus an aeeurate model is 
needed. A eomplete state-state model is an invaluable tool for testing new algorithms before imple¬ 
menting them on the aetual hardware. In order to have the model developed, a deep understanding 
of the hardware is needed and therefore the deseription of the Preeision Pointing Hexapod hardware 
is more detailed than the understanding needed to simply operate it. 

1. Actuators 

This platform is well suited for position and vibration eontrol. The aetuators, manufaetured 
by Moltran Industries, Ine. (model AFX 70N), have a stroke of more than ±5mm and ean provide 
more than 2.5^ of tilting and 10" of twisting. The aetuators ean deliver up to 40A of dynamie 
foree and up to 70A of statie foree, being adequate to eontrol large levels of vibration. The most 
important eharaeteristies, from Moltran’s datasheet, are shown on Table 2.1 and a pieture of one of 
the aetuators, as mounted, ean be seen in Eigure 2.7. 


15 





Value 

Unite 

Model 

Axial Force Transducer AFX VON 

- 

Force af dOWoc, coil af 25'’C 

100 

N 

Force af 40Wrm.s,sine, coil af 25"C 

70 

N 

Force af 40 Wdc? coil af 140^C 

70 

N 

Force af 40Wrms,sme, coil af 140^C 

50 

N 

Gain af fhe Origin 

17.5 

N 

A 

Gain af ±5mm 

15 

N 

A 


Table 2.1. Actuator Specifications. 


The actuators are powered by a custom power supply provided by CSA Engineering, Inc, 
composed by six independent switching inverting amplifiers (swifching frequency of 29KHz). The 
amplifier has negafive gain (—1) and accepfs a volfage in fhe range ±4T on fhe inpufs fhrough 
coaxial cables. The oufpufs of fhe individual amplifiers are direcfly connecfed fo fhe acfuafors. 

2. Accelerometers 

The Precision Poinfing Hexapod is equipped wifh six acceleromefers {Kistler, 8304B2 
K-Beam) mounfed inline wifh fhe sfruf’s axis, as seen in Figure 2.1. These acceleromefers have 
a range of ±2g and are linear over fhe range 0 fo 200//z. The main specifications are shown in 
Table 2.2. 


3. Disturbance Generator 


A vibrafion source is provided in order fo simulafe periodic disfurbances. This is accom¬ 
plished by using an Aura Bass Shaker AST-2B-4 mounfed on a custom adapfer. The shaker charac- 
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Value Unit 

Range ±2 g 

Sensitivity 1000 mV /g 

Frequeney Response ±5% 300 Hz 

Resolution 0.1 mg^is 

Table 2.2. Aeeelerometer Speeifieations. 

teristies are shown in Table 2.3. The shaker adapter, shown in Figure 2.8, allows vibrations along 
the z-axis, x-y axis (as defined in Figure 2.1), twist, tip and eombinations of the previous. The shaker 
ean also be mounted on the top plate, produeing vibration along the z-axis. 

The signal that drives the shaker is produeed by the DSP eontroller and sent to an AC power 
amplifier Sony STR-D Amplifier. The oufpuf of fhe amplifier is eonneefed fo fhe shaker. 

4. Floating Table and Environment Isolation 

The whole platform is sofl-mounled on a Newporf Sfabilizer^^ laminar flow isolafor, pre- 
venfing fhe ground vibrafion from inferfering wifh fhe experimenf. In order fo maximize fhe shaker’s 



Figure 2.8. Base Shaker Adapfer. 
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Value 

Unit 

Model 

AST-2B-4 


Magnet Type 

eeramie 


Power Rating 

50 

W 

Foree, Nominal, at Resonanee 

132 

N 

Weight 

1.125 

Kg 

Resonanee Frequeney 

42 

Hz 

Frequeney Range 

20-100 

Hz 


Table 2.3. Shaker Speeifieations. 


effieieney, the lower plate is mounted on the NewPort table using very soft rubber pads. Sinee these 
pads flatten due to the hexapod’s weight, they are replaeed before eaeh set of experiments to ensure 
the highest degree of isolation. 

5. Electronic Support 

Both the power supply and the anti-aliasing filters are implemented in the “Control Box.” 
This box aeeumulates the funetions of power supply for sensors, shook deteetor, emergeney shut¬ 
down, signal eonditioning and anti-aliasing filter and noise suppressor. Figure 2.9 shows the main 
funetions implemented in the Control Box. 

Power is provided by a eustom off-the-shelf power supply that provides ±5F and ±12F, 
whieh are the voltages that both the box’s eleetronies and the sensors need. 

The aeeelerometers are eonneeted by a standard DB25 eonneetor that provides them with 
-I-12F and ground. It also reeeives the aeeelerometer’s signals (« W/g), whieh are sent to the 
anti-aliasing filters. 
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Figure 2.9. Control Box Main Functions. 

A sample from the signal from the accelerometer #1 is sent to a window comparator. The 
normal operation of the hexapod should not exceed ±2g, which is the range of the accelerometers. If 
the top plate of the hexapod exceeds ±2g, the accelerometer’s output voltage will be either smaller 
than 0.5F or larger than 4.5F and thus a shock is considered to have occurred. In response to this 
condition, a dedicated circuit will disconnect the actuator and disturbance amplifiers from the 1 lOV 
rail, cutting their power supply. The power can only be restored by pressing the reset button. The 
same effect can be obtained by pressing the panic button. This mechanism is essential to preserve 
the hexapod integrity if the controller becomes unstable, especially at lower frequencies. 

The anti-aliasing filters are implemented as fourth-order switched-capacitor Butterworth 
filters. The corner frequency can be adjusted by changing the frequency of the oscillator by means 
of a potentiometer. The corner frequency is defined as fc = fosc/50. Since fhis is acfually a dig- 
ifal fiber (alfhough if is implemenfed wifh analog componenfs), one musf ensure fhaf fhere are no 
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frequencies above /c/2 when selecting fosc- The cutoff frequency was selected at 50Hz due to the 
strong nonlinearities. The filtered measurement signals are then sent to a connection panel where 
lOuF capacitors are used to reduce the noise even further. Those capacitors are needed to prevent 
the environment noise from corrupting the measurements. 

6. Controller and Host Computer 

The hexapod is controlled by a dedicated PowerPC-based board {dSpace DS1103), which 
is installed inside the host PC. A connection panel is attached to the board and provides all the 
external communication with the board. The relevant interactions are shown in Figure 2.10. 

The host computer is a PC Dell Dimension XPS T500. The controller is implemented by 
using dSpace and Mathworks’ Real Time Workshop. All filtering, disturbance generation and I/O 
operations are implemented with Simulink block diagrams. The actual control laws are implemented 
in either C or Simulink. All the data logging is performed using dSpace’s Control Desk, which al¬ 
lows the desired variables to be streamed to disk and then converted to Matlab’s proprietary format. 
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Figure 2.10. Precision Pointing Hexapod Experiment Setup. 








III. MATHEMATICAL MODEL 


A. OVERVIEW 

This chapter presents the development of a mathematical model for the Precision Pointing 
Hexapod suitable for vibration simulation and control. This model enables testing of different con¬ 
trol laws without having to implement them in the actual hardware, and thus enables the engineer 
to evaluate the actual influence of each parameter on the controller performance. 

• Section C lists the symbols used in the derivation. 

• Section D derives the equations used by the model starting from the momentum conservation 
equations, according to the Newtonian formulation. 

• Section E describes the Euler angles and Euler Parameters, explaining why the Euler angles 
were selected for the hexapod’s model. 

• Section E modifies the free-body equations to take into account the fact that the forces acting 
on the top plate are assumed to act only on the joints. 

• Section G incorporates the actual actuator’s model. 

• Section H provides the measurement equations and thus completes the model. 

B. INTRODUCTION 

Controlling a system requires knowledge of its behavior. A control system may even operate 
without a priori knowledge of the system, but it must build this information over time to perform 
successfully, as in the adaptive case. Even for adaptive controllers, a good model of the plant is 
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important in order to provide adequate parametrization. Another reason for a mathematieal model 


is given by the nature of the measurements, sinee they are usually not the quantities direetly involved 
in the requirements. In faet, the desired measurements for the platform are roll, piteh, yaw and the 
angular and linear aeeelerations. None of these quantities are measured direetly by the sensors: they 
are related to the sensors by a mapping, whieh ean be statie or dynamie. 

For the partieular problem of the hexapod subjeet of this researeh, linear and angular ae¬ 
eelerations ean be statieally derived from the aeeeleration of eaeh joint of the moving plate as a 
linear transformation. However, angular positions (|), 9 and \|/ (roll, piteh and yaw) are related to the 
aeeelerations by differential equations. To provide estimates of the desired variables, a model-based 
teehnique, sueh as the Kalman Filter, has to be implemented. The performanee and reliability of its 
results are as good as the reliability of the mathematieal model of the system. 

The rest of this seetion deseribes in detail all the steps involved in the model eonstruetion 
and presents the results. A detailed deseription of how the equations were derived is very important 
so that the same methodology ean be applied to derive models of other hexapods. 

C. LIST OF SYMBOLS, NOMENCLATURE AND DIMENSIONS 

The Preeision Pointing Hexapod was presented in Seetion II.C and only relevant details 
are ineluded in this seetion. The basie geometry of the platform is shown in Figures 3.1 and 3.2. A 
hexapod is geometrieally defined by its 12 joints (six upper joints and six lower joints). The origin is 
defined as fhe geomefrie eenfer of fhe six upper joinfs when fhe plafform is af resf wifh fhe aefuafors 
af mid-eourse (Figure 3.1). The axes of fhe eoordinafed sysfem are also depleted in Figure 3.1. 

The physieal quanfifies used fo define fhe sysfem in fhe formulalion are shown in Table 3.1. 
The definilion of fhe symbols used in fhe derivafion are in Table 3.2. 
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Figure 3.2. Hexapod—Top View. 
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Symbol Value 


.25 



r 0.0277 

0.0 

0.0 

0.0 

0.0277 

0.0 

L 0.0 

0.0 

0.0553 

100271 




Kg 

Kgm^ 


0.1271 

104271 


Description 


Radius of the circle that contains all 
the upper joints (see Fig. 3.2) 

Radius of the circle that contains all 
the lower joints (see Fig. 3.2) 

Angular separation of the clos¬ 
est lower joints, measured from 
the center of the lower joints (see 
Fig. 3.2) 

Angular separation of the clos¬ 
est upper joints, measured from 
the center of the upper joints (see 
Fig. 3.2) 

Distance between upper and lower 
planes, when the platform is at 
“zero” (see Fig. 3.1) 

Mass of the top plate 

The moment of inertia of the plat¬ 
form 

Natural frequency of the actuator 


[none] Damping factor of the actuator 

[none] Gain factor for the input to force re¬ 
lation in the actuator 

rad/s First pole of the accelerometer 

rad/s Second (and last) pole of the ac¬ 
celerometer 
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Table 3.2. Symbols Used in the Mathematieal Model. 
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D. FREE-BODY MODEL 


In order to make the problem more traetable, it is assumed that all forees and torques are 
measured with respeet to an inertial referenees system (defined in Seetion C). The most basie 
equations for a free-body model are the equations of momentum eonservation (Equations (3.1) 
and (3.2)). These qualities do not ehange if no external torque or foree is applied to the body. The 
linear (Q) and angular (H) momentum are defined as: 


Q = m\ 


(3.1) 


H = Jw, 


(3.2) 


where m is fhe mass, J fhe inerfia mafrix, w fhe angular veloeify and v fhe linear veloeify (all wifh 
respeef fo fhe fop plafe). 

The foree F and forque T aefing on fhe fop plafe are 

F = Q = m\ + m\ (3.3) 

T = H = jw + Jw. (3.4) 


Equafions (3.3) and (3.4) ean be simplified if fhe mass and fhe inerfia mafrix ean be eon- 
sidered eonsfanf, as if is fhe ease when fhe equipmenf mounfed on fhe fop plafe ean be eonsidered a 
rigid body. Wifh fhis assumpfion rh = 0 and j = 0 and Equafions (3.3) and (3.4) ean be simplified 
fo 


my 

(3.5) 

Jw. 

(3.6) 
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E. EULER ANGLES 


There are two main approaehes that ean be used to define the angular veloeity: Euler angles 
and Euler parameters (also known as quaternions). 

The first approaeh is adequate for angles within bounds, sinee this approaeh has singularity 
problems. Eurthermore, it uses transeendental funetions (sin, eos and tan), making its eomputation 
expensive. Using the azimuth-elevation-roll sequenee and aeeording to [28], the derivatives of the 
Euler angles are 




1 

tan 9 sin (|) 

tan 9 eos (|) 

0 

= 

0 

eos 9 

— sin 9 

¥ 


0 

sin(l) 

COS(|) 


COS0 

COS0 


w = M(a) w 


(3.7) 


and the singularity is at 9 = 7t/2. The platform movements are small, with all angles being less than 
15°, and the singularity will never be reaehed. Without the singularity, the only drawbaek of the 
Euler angles approaeh is its eomputational requirements. 

The other possibility is to use Euler parameters (quaternions). Aeeording to [28], they ean 
be written as 


e'l 


Wxei + Wj.ea + Wze4 

e'l 

1 

-Wxei - - Wye4 

e's 

“ ~2 

-Wyei -f w^e2 - Wxe4 

e'A- 


-Wj-ei - Wye2 - WxC-i 


with the eonstraining equation 


^1+^2 + ^3 + ^4 = 1’ 


(3.9) 


where e\, e 2 , ej, and are the quaternions. 
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This approach has two main advantages over using Euler angles: there is no singularity and 
the updates require only standard multiplications. The Euler angles approach was selected for this 
model. 

Assuming that the body and inertial frames have the same origin, a vector in the inertial 
frame °p can be transformed into the body frame using the transformation [28] 

c9c(|) —5^9 

= — celts'll/+ 5(|)5'9 c\|/ c(|)c\|/+ i'(|)i'95\|/ 5 (|)c9 °P (3.10) 

+ c(|)i'9c\|/ — i'(|)c\|/ + c(|)595\|/ c(|)c9 

= '’R/P, (3.11) 


were c stands for cos and 5 for sin. The transformation matrix is orthonormal and invertible for 
the small angles that the Precision Pointing Hexapod can produce. The inverse transformation is 

Using the Euler angles and the assumption of constant mass and inertia matrix, a simple 


model can be written as 





V 

m 

tan 9 sin (|) tan 9 cos (|) 


cos 9 

— sin 9 

silK]) 

COS(|) 

COS0 

COS0 


J“^T 



(3.12) 
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where the foree F and the torque T aet on the eenter of mass of the top plate. 


This model deseribes a body floating in free spaee and has arbitrary foree and torque as 
inputs. Therefore, this system is eontrollable as is its linearized version. The main problem with 
this model is, however, that the foree and torque are not eontrolled direetly: they are generated by 
the six aetuators. 

F. FORCES ACTING ON THE TOP PLATE 

Disregarding the gravitational effeets, there are only six external forees and no external 
torque aeting on the top plate (through the joints). A mathematieal relationship must be found that 
ean translate these forees into the resultant foree and torque. It is assumed that the aetuators are 
mass-less and the forees generated by the aetuators are parallel to the aetuator’s axis and aet only at 
the joints. 

The foree aeting on the eenter of mass is the summation of all forees aeting on the body. 
Defining L, as the foree aeting on the joint i and e = [ei as the foree produeed by the 

aetuators, the resultant foree is 

F = 


where e, is the foree produeed by the aetuator i, b, is the position of upper joint i and B, is the 
position of the lower joint i. 


6 

IC 

i=l 

^ b,-B, 


(3.13) 


b, = p + r, 


= p + *R„ro,,- 


(3.14) 
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Defining iT, = the foree F ean be written as 

F = 

i=i 

I 


7x 7i 53 54 55 56 


= 5(p,a)e. 


(3.15) 


The torque ean be evaluated in a similar manner. Eaeh foree has a assoeiated torque when 
translated to the eenter of mass and these torques ean then be summed to produee the total torque 
T. First, define a veetor r,- as 


r, 


Rr,r„ 


(3.16) 


where p is the position of the eenter of mass of the top plate and b, is the position of the joint i. 


Therefore, the overall torque ean be written as 


T 



b, -B, \ 

b, -B, \ 

|b,'-B;||y^ 


'F(p,a) e. 


(3.17) 


where 5, = r,- x 


b,-Bi 

[bi^ 
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The results of Equations (3.15) and (3.17) provide a linear relationship between the input 


forees and the resulting foree and torque. By using these results, the plant model ean be rewritten as 

P 

V 

d 

w 

where the parameters of iT (p,OC) and ‘T (p,(x) were dropped for eompaetness and are implieit. 

G. ACTUATORS 

The foree of the aetuators eannot be eontrolled direetly. The eontroller’s output is a num¬ 
ber that is eonverted to a voltage/eurrent and then sent to the aetuator. The aetual foree produeed 
depends on the aetuator’s internal dynamies. 

The aetuators installed on the Preeision Pointing Hexapod are voiee-eoils, whieh are ba- 
sieally a eoil in the inside and a magnet on the outside (Figure 3.3). A eurrent flowing in the eoil 
generates a magnetie field that interaets with the magnetie element ereating the foree, whieh is linear 
for small displaeements. Both a spring foree and a viseous foree must also be taken into aeeount. 
Therefore, the equation relating the foree generated by the aetuator to the eurrent flowing through it 
and the position is 

F{t) = Kii{t) — Kxx{t) — bx{t) — max{t), (3.19) 

where Ki is the gain from eurrent to foree and b, nia and Kx are the viseous frietion, moving mass 
and spring eonstant of the aetuator. 


m'^ 

M(a) w 

J-^Te 


(3.18) 
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Figure 3.3. Actuator Diagram. 


The force produced by the actuator must be written in terms of the 
tern. The displacement .r is the distance traveled by the joint along the unity 

Xi = J'i 

= (p + r,-bo,;)^5; 

= (^p + ^Roro,, -bo,;) J,-. 

Since b^,; = ro,;, this becomes 

^ = (p + "Robo,-bo,,y5,-. (3.20) 

The velocity v, of the joint i is the sum of the velocity of the center of mass and the velocity 
of the joint with respect to the center of mass. Since the latter is produced by the angular velocity 
of the top plate, the resulting expression is 

Vi = (v + wxr^j^jTi 

= (\ + yfx‘’R„ro,t^ Ji. (3.21) 


inertial coordinate sys- 
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Similarly, the acceleration of the joint i is the acceleration of the center of mass summed to 


the centripetal acceleration: 

a,- = a + w X (w X r,-). (3.22) 

Projecting a, it on the struts’ axes leads to: 

Ui = (a + wx (wxr,))^jr,- 

Ui = (^a + wx ^wx^R^ro,;^^ iT,-. (3.23) 

Substituting x, x and x into Equation (3.19), the force on joint i is 

Fi{t) = Kiii{t) - + ^Robo,,' - bo,;^ +b(\ + y/x + 

(^a + w X X ^ ^ 5,-, (3.24) 

assuming it is parallel to the actuator’s axis. All the terms of this equation are known states except 
for the input /,(t). 

It is important to mention that the output of the controller is not a current, but a voltage. This 
voltage is sent to a DC amplifier and its output goes to the actuators. The amplifier’s bandwidfh is 
200Hz and fhe fime consfanf of fhe acfuafor’s LR circuif is Im^. Therefore, fhe model developed in 
Ibis secfion is valid only up fo 200Hz- If higher frequencies are needed, fhen fhe elecfrical dynamics 
musf also be modeled. 

1. Model With Actuators 

The force Equafion (3.24) can be used in Equafions (3.12) fo generate a more complete 
model. To do fhaf, Equafion (3.24) musf be wriffen in a more compacf form. Defining 

hi = (^Kjc + ^Robo,,' - bo,;^ +b(\ + v/x + 

(^a + w X X ^RoCo,,-^ ) ) ^'’ 
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h^[hi ...hef 


(3.26) 


and 


u= [/i ... kf. 


Equation (3.12) can be rewritten as 


p 


V 

V 


m'^ 

a 


M(a) w 

w 


r^'Te 




E 1 ...F 6 

M(a) w 




Fi...Fe 


iyWi'-h) 

M{a)w 
[j-i‘r(^,u-h) 
with h and u as defined in Equations (3.25), (3.26) and (3.27). 


(3.27) 


(3.28) 


H. ACCELEROMETERS 


All hexapods used for vibration control reviewed in the literature have some kind of sensor 
to measure the vibration built into the strut. The Precision Pointing Hexapod is not an exception 
and has accelerometers built inline with the actuator’s axis. Section ll.C.l (p. 15) describes the 
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mounting bracket and the accelerometers, while Table 3.3 shows the relevant specifications of the 
accelerometers. 



Value 

Unif 

Range 

±2 

g 

Sensifivify 

1,000 

mV/g 

Resonanf Freq. (nom.) 

1.2 

kHz 

Freq. Response (±5%) 

300 

Hz 

Finearify 

±1 

%FS 

Transverse Sensifivify (nom.) 

< 1 

% 

Phase Shiff (lOHz, nom.) 

2 

degree 

Resolution (0.5 fo lOOHz) 

0.1 



Table 3.3. Accelerometer Specifications. 


As seen in Section G, the model is only valid up to 200Hz due to restrictions of the actuator’s 
model. Since 200Hz is within the accelerometer’s bandwidth, their dynamic can be considered just 
a constant. It is also important to realize that the accelerometers mounted on the struts measure DC 
level, and this model does not limit this behavior. 

The accelerometers measure the acceleration of the joint projected on the actuator’s axis. 
The measurement can then be written as 


ai = ^ 9 i 

b; -B; 


(3.29) 


where a, is the measurement of accelerometer i and a, is the acceleration at the joint i. This acceler¬ 
ation was evaluated in Equation (3.23). 

Using the above result and defining as fhe acceleromefers’ consfanf gain, fhe measure- 
menf equations become 


m 


ai 




ae 



Ka, (a + wx (wx*Roro,i))^iri 


Ka, (a + w X (w X *Roro,6)) ^ ^6 


(3.30) 
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Therefore, the last equation is the measurement equation for the state-spaee model de- 
serihed in Equation (3.28), whieh ean he used for vibration simulation purposes. 
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IV. VIBRATION SUPPRESSION 


A. CHAPTER LAYOUT 

This chapter contains the details of the Adaptive Disturbance Canceller. First, Section B 
derives the Multiple-Error LMS method. 

Section C presents an overview of the Adaptive Disturbance Canceller. The mathematical 
details, including a convergence analysis, are presented in Section D. 

Note that the Adaptive Disturbance Canceller isolates the vibrations at the assigned fre¬ 
quencies only. As a result. Section D.3 presents a study of the behavior of the controller in the 
presence of unassigned frequencies on the error signal. 

Section E presents simulations results obtained with the Adaptive Disturbance Canceller. 
The model derived in Chapter III was used to simulate the PPH in Section E.2. The Adaptive 
Disturbance Canceller is a SISO controller and multiple controllers might interact with one another 
on a distributed controller configuration. Thus, full nonlinear model simulations are of special 
interest. 

B. MULTIPLE-ERROR LEAST MEAN SQUARED 

Advances in Digital Signal Processing techniques led to several estimation techniques and 
many of these are also valuable for noise cancellation. Techniques commonly applied in noise 
cancellation are usually based on adaptive algorithms. The general approach is shown in Eigure 4.1. 

Among adaptive filter implementations used to cancel the disturbance signal d [n] one of the 
simplest and best understood is the Eeast Mean Square method (see [29] for a detailed description 
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Figure 4.1. Noise Caneellation Problem. 


and speeifie properties). This filter assumes an all-pole plant {H) and minimizes the eost funetion 

7 = E{|eHn, (4.1) 

where E{«} is the expeetation operator. Additionally, the LMS is a first order method, whieh 
eontributes to its simplieity and redueed eomputational effort in eomparison with more sophistieated 
algorithms. 



Figure 4.2. An Adaptive Vibration Controller. 


The vibration problem has one additional eomplieation when eompared to the noise ean- 
eellation problem: the transfer funetion H is after the filter. In faet, the filter provides the inputs to 
the plant H. An updated diagram, refieeting the vibration eontrol problem is shown in Figure 4.2. 
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It is important to realize that Figure 4.1 and Figure 4.2 only differ on the position of the plant and 


the eontroller: they are switehed. Unfortunately, the very faet of switehing the plant and the filter 
prevents a straight LMS from being used for eontrol. The LMS filter has a weight update equation 
in the form 


W[u+1] = WH+qe*[u]u[u], (4.2) 

and sinee the plant is plaeed after the filter, a ehange on u [n] will not refleet on e [n], but on e [u + 1] 
(the LMS assumes an all-pole plant). Therefore, e [n] does not depend on u [n] and the LMS method 
eannot be applied direetly to the vibration eontrol problem. 

Widrow[30] and Burgess[31] solved this problem working independently in 1981. Their 
method used an estimated model of the plant and is ealled Filtered-X LMS. This algorithm inherits 
most of the eharaeteristies of the simple LMS scheme. Unfortunately, the Filtered-X LMS method 
does not address one important issue: MIMO systems. Elliot, Stothers and Nelson generalized the 
original Filtered-X LMS in 1987 [12], presenting the Multiple-Error EMS. The new algorithm was 
able to control MIMO plants, but still possesses the restriction of having a single signal correlated 
with all disturbances. In essence, the Multiple-Error EMS transforms the MIMO plant into a SIMO 
plant, with a single input {x{n]) and several outputs {eo{n]) (see derivation in Section 1). The 
structure of the controller proposed in [12] can be seen in Eigure 4.3. 

Both the Multiple-Error EMS and the Eiltered-X EMS approaches, as well as their variants, 
seem to be the most common active vibration controFisolation algorithms in use today. They are 
among the least expensive in terms of computational requirements and are comparatively simple to 
implement. The Multiple-Error EMS, due to its widespread acceptance, is used in this work as a 
benchmark and is presented next. 
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Figure 4.3. LMS for Vibration Control. 



Figure 4.4. Multiple-Error LMS. 


1. Derivation of the Multiple-Error LMS 


The Multiple-Error LMS assumes that the MIMO plant H (Eigure 4.4) ean be modeled as 
a Einite Impulse Response (EIR) filter or order M with Q outputs and P inputs. Thus, the q*^ output 
[n\ due to the input p is given by: 


M 

yq,p W] = L H “P [n-m], 

m=0 


(4.3) 
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assuming that the Q disturbances dq [n] acting on the system are correlated to the signal ;c \n], which 
must be measurable. 

Further, assuming W to be a matrix of dimensions S x P, the signal u [n] can be generated 
as 

5-1 

“p [n] = [”] (4.4) 

.s=0 

Substituting the input value derived in Equation (4.4) into Equation (4.3), the output of the 
plant at time n is 

P-l M 

yq [«] = L L ^q^p C'”] “p “ '"I 

p=0m=0 

P-l M 5-1 

= X 12 Vp H 12 ^■'^P [n-m]x[n-m-s\ 

p=0m=0 A—0 

P-l M 5-1 

= 12 12 Y^Kp[p\ws,p['^-m\x[n-m-s]. (4.5) 

p=0m=0.s—0 

The error can be then computed as 


Cq [n] = dq [n] - Jq [n] 

p-l M 5-1 

= <5?^ [«] - XI 12 12 ^q^P w ^*:P [n-m\x[n-m-s\. (4.6) 

p=0m=0.9=0 

Eollowing the same reasoning that is used in the derivation of the EMS method, W must be 
updated minimizing the cost function 


7 = E 


e-i 

L 4 




(4.7) 
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By using the steepest deseent and as the learning rate, the weight ean be adjusted as 


follows: 


Ws,pW\ = Ws,p[n]-n 


^Ws,p [n] 


= Ws,p[n\-fi- 




dws,p [n] 


Q^dE{el[n\} 

= w,^p[n\ -r-p. 

to 


Unfortunately, eomputing the expeetation E [n ]} in real time is impossible, but the same 
approximation used in the derivation of the standard LMS algorithm ean be made here: 


E{el[n]}^el [n]. 


By using the above approximation and proeeeding with the ealeulation of Ws^p [n + 1]: 


w,,p [n + 1] « Ws,p [«] - f/ 2, ^ r 

q=0 

^ r 1 o ^ r 1 
Pi w,^p[n]-2iil^eq[n]^^ 

q=Q OWs,p [« 


(4.10) 


Next, substituting the error Cq [n] from Equation (4.6) into Equation (4.10), Wg^p [n-\- 1] ean 


be evaluated as 


Ws,p [n+l] ~ Ws,p [n] + 




1 

. A n\ __ L _^_ 


(4.11) 
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Assuming that the weights ehange slowly (ju small), then, Ws^p [n — m] ^ Ws^p [n] and the 


previous equation ean be simplified to 


Ws,p [n+l] « Ws,p [n] + 


hLj /[«—mkf/t—m— 

r 1 y“l}'=0^m=0^s'=0 q,p''- ^ J L ij 

Ws,p [n] + 


nLj /[oljrfn 

. ^Q-l r 1 V /^=0 5^=0 ^^P '- ^ s',p''-^ L 

2a'i:jLo «d«] ^- 


Ws,p [n] + 


o y'G"! r 1 f v/’—1 v'M } r 1 

^q[^\ ^-=0 ^y=0 Kp' [H 


x\n—m—s 


(4.12) 


Note that the derivative 


tion (4.12) eollapses to 


9«'y.p'H ; 


'dws[p[n\ 


is different from zero only if p = p' and s = s'. Therefore, Equa- 


e-i M 


Ws,p [u + 1] Ws,p [n] + 2q ^ [n] ^ hg^p [m]x[n-m-s]. 


(4.13) 


Now, defining rq^p [n] as 


'q,p 


fhen, Equafion (4.13) beeomes 


^=0 m=0 


M 

[n] = ^ hq^p[m]x[n-m] 
m=0 


e-i 


(4.14) 


Ws,p [u + 1] ~ Ws^p [n] +2pY^eq [n] Vq^p [n-s\. 

q=0 


(4.15) 


Equafion (4.15) is fhen used fo define fhe weighf updafe equafion for fhe Multiple-Error 


EMS: 


e-i 

Ws,p [n+f]= Ws,p [n] + 2q ^ [n] rq^p [n-s], 

q=Q 


(4.16) 


One imporfanf eharaeferisfie of fhe Multiple-Error EMS mefhod is fhaf a planf model is 
needed fo generafe rq^p [n]. Aefually, fhe ferm Y,m=o^q,p [p]x[n — m — 5 ] ean be inferprefed as being 
fhe q'^ oufpuf of fhe planf when fhe p'^' inpuf is x [u — 5 ] and all fhe ofher inpufs are zero. 
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C. PROPOSED METHOD 


For a time-invariant or slowly varying tonal disturbance, having a full model of the plant 
or a controller designed to work over a frequency band is an overkill. The system is only excited 
at a particular frequency and thus, the transfer functions are reduced to single complex numbers. 
Assuming that the plant is linear and stable, it is possible to generate any arbitrary sinusoidal at the 
output of the system and thus cancel any sinusoidal disturbance. A controller using this approach 
was presented in 1998 by Bertran and Montoro [18]. The proposed controller needs, for the studied 
case, only a few assumptions: a stable linear SISO plant, tonal disturbances with known frequencies 
and the plant without zeros at the frequencies of interest. A block diagram describing the controller 
is shown in Figure 4.5. 

The controller proposed can be summarized as follows: assuming that the plant H is linear 
(see Figure 4.5), then, for any sinusoidal signal d[n] with frequency (Oc it is possible to find a 
sinusoidal input x [n] such that y [n] = —dn (if (e^®^) ^ 0). Usually one would write the input as 

x[n] = AcOs((Oc?^ + Pa:)- (4.17) 

Several algorithms which compute the optimal value of x„ by minimizing e[n\=y [n] + d [n] 
are available, but almost all assume that y [n] is a linear combination of the parameters. The input 
Equation (4.17) can be changed to the equivalent form 

x[n] = acos((Oc«)+ ^sin(c0c«)- (4.18) 

Next, assuming that 9-C (e^®‘^) = ae^P, the steady-state output y [n] can be written as 

y[n] = aacos(c0c?^ +P) +^ctsin(c0c?^ + P)- (4-19) 

Using this form, the output y[n] is linear in the parameters a and b and most adaptive 
algorithms can be used to find a and b. Berfran and Monforo used an updafe similar fo fhe one 
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Figure 4.5. Adaptive Disturbanee Caneeller. 


defined in the LMS method: 


a[n-\-\] = n[n]+/re[n]eos(cOc?^) (4.20) 

b[n-\-\\ = b[n\-\-jJLe[n\?,m{G)cn), (4.21) 

with 

e[n]=y [n] + d[n]. (4.22) 

D. MATHEMATICAL ANALYSIS 

1. Stability 

Bertran and Montoro derived the stability analysis for the setup presented in this paper, but 
the partieular plant assumed is not appropriate for the Stewart Platform problem. Therefore, a more 
general approaeh is presented in the following diseussion. 

The disturbanee d [n] ean be written as 

d[n] =AtlF(e-^'“^)e-^'“^”, (4.23) 
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where cOc is the frequeney of the disturhanee. 

Assuming q small, A [u] ehanges slowly and the output y [n] ean be written as 

y [u] A [n] ^ (e^'“^) (4.24) 

Thus, sinee e[n\=y [n\ -\-d[n], the error is 

e\n] = d [n] + y \n\ 

= (A + A[u])l//(e-^'®-)e^'“-” 

= A[u]l/f(e-''®^)e-^'®^". (4.25) 

Using the weight update equation and the definition of A [u] = A+A [u], the update equation 
for A [n\ is given by 

A [u + 1] = A [u] + qe [n] e“-^®‘'”. 

Adding A to both sides of the previous equation leads to 

A [u + 1] = A [u] + qe \n] e“-^®^”. (4.26) 

Also note that algorithm eonvergenee may be proven if |A [u + 1] | < |A [u] |, for all integer 
values of n. Therefore, 

|A[u+ 1]|^ = A* [u+l]A[u+1] 

= (A* [n] + qe* [n] e-^'®^") (A [n] + qe [n] e”-^'®^") 

= I A [n] I V A* [n]^ie [n] e”-/'®^" + A [n] e^'®^” + \^i\^\e[nf. (4.27) 

Now, using the error value from Equation 4.25, the above expression is redueed to 

I A [u + 1] I ^ = |A [u] I ^ + A* [n] liA [n] d{ (e^®^) e^'®-"e“^'®^” + 

A [nYA* [n] di* (e^®^) + 

Iql^A* [n] d{* (e-^®^) e--^'®^"A (e-^®^) e-^'®^". 
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Applying complex conjugate properties, the final expression for |A [n+ 1] ^ is given by 


IA [n + 1] I' = \A[n]f + lRe{A [nYA* [n] (e^'“^)} + |p|" | A [n] |" | iff (e^'“^) |' 

= |AH|^(l + 2/?e{A/iff(ef“^)} + H^|^(e-''“0|^)’ (4-28) 

and it can be concluded that lim„_><x> | A [n] | ^ = 0 provided that 

1 + 2Re {n9{ (ef“^)} + \n\^ |iff (ef“^) < 1, 

or, equivalently, 

2Re (e-'®^)} + |qp |iff (ef®^) < 0. 

The weight A [n] converges to —A, as A [n] = A + A [n]. Thus, the error converges to zero. 
Assuming that |qiff (ef®^) | << 1, the stability criterion can be simplified to 

Re {qiff (e^®^)} < 0. (4.29) 

As a guideline, q can be selected by using 

iff* (ef®A 

f'=^ -(4.30) 

10|iff(ef“0r 

where the magnitude of the plant is not critical and the phase must be known within ±7t/2. This 
guideline is not strict and the model of the plant does not need to be precisely known. 


2. Optimal Weights 


Assuming that the controller converges, then y [n] = —d [n]. Assuming that the disturbance 

d[n] is 


r/[n] =D,ef(®^"+^4, 


(4.31) 
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then, the output y \n] can be evaluated as 


e[n] = y [n] + d [n] 


The convergence of the algorithm ensures that the error converges to zero and A [n] con¬ 
verges to a constant. Therefore, 


e[n\ = y [n] + d [n] 


0 = Al//(e-^®^)e^®^'’ + D, 


j&A J&ctr I p. j{iacn+yc) 


Defining d{ (e^®^) = ae^P, the previous expression becomes 


= + aAe^'(®^”+P) 


n e^'(®^”+T^-) = -(xA*e^'(®^”+P) 


lEgAYc-p+’t) 


(4.32) 


where A* is the optimal weight. 

The real coefficients a* and f?* are obtained by rewriting Equation (4.32) in the real form 

a* = -^cos(Yc-P) 


= 


■sin(Yc-p). 


(4.33) 


3. Crosstalk 


It is important to verify how the algorithm behaves in the presence of uncontrolled fre¬ 
quencies. For such an analysis the disturbance is assumed to contain two complex exponentials at 
frequencies 00^ and (0„c- Thus 


d [n] = 


(4.34) 
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where (Oc is the disturbanee to be suppressed and co„c is not to be suppressed. 


It is assumed that € N (where N denotes the set of Natural numbers) sueh that A [n] = 
A[n + N] and 


< ytg 4 e and 


The input of the plant can be written in the complex exponential form as 


x[n] =A[n]e^'“^". 


(4.35) 


Taking the DFT of x [n] and applying the frequency shift property, X (k) becomes 


X{k) = AL{k-kc). 


The output y \n] is 


y\n] = X\n] *h\n]. 


and therefore, {y [«]} is 


‘y{k) = 

= Al{k-kc)9{(&^^^^y 


(4.36) 


The weights A [n\ are updated as 


A[n+1] = A[n]+qe[n]e 


(4.37) 
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Since e{n\=y [n] -\-d{n], Equation (4.37) can be rewritten as 

A [n + 1] = A[n\-\- jji{y[n]-\- d[n]) 

= A[n]+ lue-j^^’^y[n] + qe--^'®^” + D,e-^'®-"+^'^-) 

= A[n] + lue-j^^’^y [n] + qD,e^'^^e-''®^'’e-^'®^” + qD„ce-/'^-e^'®-"e--''®^" 

= A[n] + fue~-'^y[n]+ fjDce-'"^^^ + fjDncQ-'^e-^^’"'^. (4.38) 

Applying the DFT to Equation 4.38, leads to 

®)F‘r{A(n+1)} = ®)r‘r{A[n]}+q!D)r‘r {e“-^'®^”y[n]} + 

©(FT {1} + ©(FT | e-^'^“” |. 

Using DET properties, the previous equation becomes 

• 27CF 

Jl{k)eJ^ = AL{k)+ij‘X{k + kc)+ijDceJ^^5{k)+ijDnce^'^'^5{k-ks). (4.39) 

By substituting “X (k) from Equation (4.36) into Equation (4.39) leads to 
AL (k) = J4 (F) + (F) di ^ 5 (F) + 5 (F — F 5 ) . (4.40) 

The final expression can then be obtained by grouping the terms in AL (F) 

J^{k) (ej^ -I-^d{ ^e^'(T+®c)^^ = qDce^'^^5(F) + 

ljD„,&^^'-5{k-k8). (4.41) 

Equation (4.41) shows that the weights have components only at F = 0 and k = k^ and 
thus the error only has components at (0^ and (0„c- Therefore, the controller does not create new 
frequencies. Influence of the controller in the error at co = co„c is evaluated by replacing F = F 5 in 
Equation (4.41) which leads to 

Al (F 5 ) ^e^“^ — 1 — fjdf ^e^( " ^ ^ = qDce^^'"5 (F 5 ) + (F 5 — F 5 ) 

= (4.42) 
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Thus, 


mh) = 




jjnt 


qJ-tt — 1 — fj^ (e-/“"£') 


Similar results may be obtained by evaluating Equation (4.41) at = 0 


(4.43) 


(0) (e^'^ - 1 - qi/C ^e^'( ^+“^)) ) = qD^e^'^'^S (0) + qD„ce-^'^"^5 (0 - h) 

= (4.44) 


Thus, 


MO) 




(4.45) 


Using the fact that ‘X{k) = J^{k- kc) ^ ‘E{k) = ‘y {k) + ‘D{k) can be evaluated 


as 


£ (fc) = ‘y{k)+ 5 (^ - ^s) + 5 {k) 

= J^{k — kc) ^ + Dnce^^'^""?>{k — knc) + DcQ^’^'h{k — kc) 


( D 

I — J—. — -b{k — kc) + 


DncQ^'^''‘^^{k - knc) + DcM‘^5 {k - kc ). 


• 27 C*s 

qJ^T — 1 — fjfkC M^"") 


+ 


(4.46) 


At this point. Equation (4.46) shows that the error has components at COc and (One only. 
Eurthermore, the disturbance at (One does not affect the error at (0^- Accordingly, the disturbance 
at cOc does not affect the error at the frequency (One- Thus, the controller does not generate new 
frequencies. 

Evaluating Equation (4.46) atk = knc the error due to the noncontrolled disturbance is 

•E{knc) = .2.,, -i;f(et®-)+D„,e^'^-. (4.47) 

qJ n _ 1 _^j;/'(ei®nc) 
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Assuming that So is small {k^ « N), then e " 1 and Equation (4.47) can be simplified 


to 


£ {knc) = 




Jjnc 


iKk^ 


qJ N — 1 — (ef®nf ) 


(e-^®”) + Dnce-^'^"‘ 


—fjiH (e-t®"f) 
( 1 - 


^ +D„ce-^^“ 


1 






0 . 


(4.48) 


Equation (4.48) shows that the error tends to zero when the assigned frequency tends to 
the exact disturbance frequency (6(o tends to zero). This important result shows that there is not 
discontinuity in the solution in the vicinity of the true solution, provided (0„c is not a zero or pole of 
the plant. 

However, assuming 6 q) is sufficiently large and (e^®"^) | << 1, leads to 


iiD 

<r(k ) - ^ _ 

^yf^nc) — ,2Ttk^ 

qJ^t — 1 — (ef®«f) 




fjDfic^ 


Jlnc 


e-' " — 1 


(e-^®"^) + 


1 + ^ (e^®-) . 


(4.49) 


e-^ « — 1 

So, using \fjfH (e-^®"^) | << 1 and assuming that ||e^^“ — 1|| 1, Equation (4.49) may be 

approximated by 


£ {knc) ~ Dnce-’'^"‘' 1 + 


pi p,jync 

i^nr^ 


.iKk^ 


H (e-^®"") 


(4.50) 


The above result shows that the controller has very little impact on the component at (0„ 
when this frequency is far from 00^ and q is small. 


54 



E. SIMULATIONS 


Simulations are conducted to verify the mathematical development and to verify whether 
the assumptions made are reasonable. Recall that stability conditions were derived assuming a 
single frequency on the error signal, and thus it is important to verify the behavior of the controller 
if other frequencies are present in e [n]. Additionally, the Adaptive Disturbance Canceller is a SISO 
controller; therefore, it is necessary to verify if it can control a hexapod, as it is a strongly coupled 
MIMO plant. 

1. Single-Input/Single-Output System 

The first simulation is intended to evaluate the algorithm convergence when the distur¬ 
bance consisted of a single tone. The method was implemented in Simulink as described by Equa¬ 
tions (4.20) and (4.21). The plant used was a second order system with a natural frequency of \0QHz 
and damping of 3% and its output was sampled 1,000 times per second. The controller also ran at a 
\KHz sampling rate. White noise was added to the error signal, with SNR = 20dB. The results are 
shown in Fig. 4.6. 

Results show that the controller works as expected, even for such a high noise level. In 
addition, note that the controller does not affect the system at uncontrolled frequencies. 

The next simulation studies the controller behavior when two tones are present in the dis¬ 
turbance signal (fc = 33Hz and /„c = AlHz), where 33Hz is the disturbance to be suppressed. Fig¬ 
ure 4.7 shows that the controller cancels the assigned frequency and does not change the unassigned 
one {AlHz)- The weights exhibits components at / = QHz (from fc), f = lAHz and / = %QHz- The 
latter two values are predicted by theoretical results, since |/^ — fc\ = 14//z and fc + fd = ^OHz- 
Figure 4.8 shows the performance of the algorithm when noise is added to the sensors {e [n]). 

Theoretical results predicted that this controller cannot be applied directly when the phase 
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Figure 4.6. Single Noisy Tone. 

of the plant is | + Ikn because 11 + /rae-'P | = 11 + jpa sin( p) | > 1. In order to verify how sensitive to 
this parameter the method is, the plant was changed so that its phase at fc is almost 90°. Figures 4.9 
and 4.10 show the resulting performance of the controller. As expected, the maximum |^| is very 
small and thus the convergence rate is very slow. 

The last SISO simulation tests the performance of the algorithm when several controlled and 
some uncontrolled frequencies are present. Figure 4.11 demonstrates the results for four controlled 
(10, 14, 34, and 36Hz) and three uncontrolled (13, 35 and 40Hz) tones. Noise was added to the error 
signal. It is relevant to mention that among the controlled frequencies there were 34, 35 and 36Hz, 
where only 34 and 36Hz are to be suppressed. The goal was to verify how selective the algorithm 
is. Results show the tone at 35Hz is left untouched while the tones at 34 and 36Hz are attenuated. 
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Figure 4.7. Case B, No Noise. 
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Figure 4.8. Case B, Noisy. 
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Figure 4.9. Case C, No Noise. 
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Figure 4.10. Case C, Noisy. 
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Figure 4.11. Case D: Four Controlled and Three Uncontrolled Tones. 
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2. 


Multiple-Input/Multiple-Output System: The PPH 


The Precision Pointing Hexapod is a highly coupled system as it is impossible to move any 
node separately from the others. The Adaptive Disturbance Canceller was developed as a SISO 
controller and either the controller, must be generalized to a MIMO approach, or a multiple-SISO 
approach must be taken (also called decentralized control). Previous attempts in using decentralized 
control had mixed success. Beavers [11] did not obtain a good performance for higher frequencies 
and the performance for lower frequencies was poorer than the one obtained by using the MIMO 
methods (Multiple-Error LMS, Repetitive Control and Clear Box [4, 14, 13]). Good results using 
distributed control (six independent SISO controllers) were obtained with some controllers [27, 32]. 
Since previous work exhibited mixed results, the mathematical model developed in Chapter III 
(p. 23) is followed to verify if the coupling is a limiting factor for the use of the Adaptive Disturbance 
Canceller on hexapods. 

The simulation setup was exactly the same as that used in Section E (p. 55), but a multiple- 
SISO implementation was used and the plant was changed to an S-Function implementing the model 
of the hexapod according to Equations (3.28) and (3.30), on pg. 36. Simulations were made to excite 
the piston (bouncing), tilt/tip, twist and shear modes (see Eigure 4.12). 

The first case simulated shows the performance of the controller for the vibration axis co¬ 
inciding with the z-axis and no rotational component (piston mode), with a frequency of 40Hz- The 
results are shown in Eigures 4.13 and 4.14. As one can see, the decentralized controller performed 
reasonably well and no competition among the controllers was seen in steady-state. The nonlin¬ 
earities in the plant model are responsible for the harmonics seen in the error signal observed in 
Eigure 4.13(a). 
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Figure 4.12. Hexapod modes [5]. 


The seeond ease simulated was the twist mode in whieh the aetuators needed to move 
symmetrieally in order to eaneel the vibration. As Figures 4.15 and 4.16 show, the performanee was 
as expeeted and was the same as that obtained in the piston mode. 

In the next ease, a more general vibration disturbanee was simulated by making eaeh eom- 
ponent of the signal dn have a different amplitude (namely 1, 3, 2, —4, 3 and —2m/s^). Figures 4.17 
and 4.18 show the eontroller performed as well as in the previous two eases. 
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The last simulation added a seeond tonal {50Hz) disturbanee to the previous ease (ampli¬ 
tudes 2,1,6,—12,1 and —5m/s^). This disturbanee was not assigned to the eontroller and therefore, 
was not being eontrolled. The amplitudes were the same as those for Case 3. It ean be seen in 
Figure 4.20 that the eontroller performed as expeeted. 
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Figure 4.14. PPH Simulation, Case 1, Noisy {SNRd = 50dB) (Accelerometer#!). 
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Figure 4.15. PPH Simulation, Case 2, No Noise (Accelerometer #1). 
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Figure 4.16. PPH Simulation, Case 2, Noisy {SNRd = 50dB) (Accelerometer #1). 
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Figure 4.17. PPH Simulation, Case 3, No Noise (Accelerometer #1). 
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Figure 4.18. PPH Simulation, Case 3, Noisy {SNRd = 50dB) (Accelerometer #1). 
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F. COMMENTS 


The simulations verified the theoretieal results and the Adaptive Disturbanee Caneeller ex¬ 
hibited good performanee even when used in a distributed eonfiguration and in the presenee of 
strong nonlinearities. It is also important to mention that the vibration suppression was not de¬ 
graded by the presenee of noise in the measurements. This result was eonsistent for both the SISO 
and the MIMO simulations. 

The simulations also eontirmed the theoretieal results about the unassigned frequeneies: the 
eontroller does not affeet the error at frequeneies different from the assigned ones, indieating that 
this eontroller is very seleetive. 

It is important to reeognize that the properties mentioned in the last two paragraphs are 
interrelated: one exists beeause of the other and viee-versa. Sinee the eontroller uses a synthetie 
signal as input, it is very insensitive to noise. But a synthetie referenee signal assigned to an ineorreet 
frequeney will prevent the eontroller from eaneeling the disturbanee. 

As expeeted, the eontroller did not generate any new frequeney in the SISO simulations, 
whieh used a linear plant. The Preeision Pointing Hexapod model used in the MIMO simulations 
eonsiders the geometrie nonlinearities of the hexapod. These nonlinearities generated harmonies 
that appear over the entire range of frequeneies, although the harmonies are very small beeause 
the displaeement of the top plate is small and thus the angles are small. Nonetheless, this effeet 
is expeeted to be amplified if other nonlinearities sourees are present, sueh as statie frietion and 
baekslash. 

Note that one important limitation of this method is the requirement for the preeise fre¬ 
queney of the disturbanee or the performanee will suffer eonsiderably. This effeet was already 
observed by Tsuei et al. in [33], whieh presented a modified version of the Filtered-X LMS with 
synthetie input. 
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V. COMPUTATIONAL COST ESTIMATION 


A. OVERVIEW 

This chapter compares the computational cost required to implement the Adaptive Distur¬ 
bance Canceller and the Multiple-Error LMS method. Both computational complexity and memory 
requirements are compared. 

Section B defines the methodology followed to estimate the computational cost. Section C 
and Section D present the computational cost for the Multiple-Error EMS and the Adaptive Distur¬ 
bance Canceller respectively. 

Section E compares the computational cost and memory requirements for both controllers. 
One example is presented in which a hexapod is used to suppress three tones with one harmonic 
each. 

B. BASIC CONCEPTS 

Vibration control requires large amounts of computational power. The frequencies involved 
are usually higher than the traditional control methods; therefore, the sampling rates are also higher. 
One of the main purposes of this research is to study a vibration method that requires minimal 
computational resources. The controller, after being implemented, should consume the least amount 
of electrical power possible. Additionally, the computer system should require as little space as 
possible, since electrical power, weight and space are very expensive in space applications. Ideally, 
a single low-power microcontroller, with no external memory and no external A/D converters should 
be used. 
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Floating point operations used to be the most expensive operations; eonsequently, the eom- 


plexity of an algorithm is traditionally measured by eounting floating-point operations. It was not 
uncommon for a simple floating-point multiplication to be emulated in software or to need a ded¬ 
icated external chip to be executed. The advances in microelectronics changed the picture: CPUs 
have evolved. The new DSP CPUs can perform one floating-point multiplication as fast as they can 
perform one floating-point summation or even integer math. Table 5.1 shows the specifications of 
an embeddable microcontroller/DSP (Hitachi’s SH7410, a CPU for embedded devices). This mi¬ 
crocontroller CPU has a DSP unit embedded, interrupt controllers, memory, A/D converters, timers, 
etc. The power consumption is very small (0.5W versus about 60W for an Intel’s Pentium III or 
an AMD’s Athlon). Therefore, if a vibration control solution can fit in a few dozen kilobytes of 
ROM and does not need more than a few kilobytes of RAM, the savings in space, weight and power 
are very significant. As an even simpler alternative, the controller could be converted to use integer 
arithmetic or floating-point emulation and use a CPU that does not support floating-point. These 
are more easily available and come with a wider array of features (some featuring eight 12bits A/D 
conversors and eight D/A conversors), making it possible to implement the controller with very few 
components and little power consumption. 


Voltage 

3.3V 

Operating current 

110 to 210mA 

Operating power consumption 

0.36 to 0.691T 

Standby current 

5/rA 

Standby power 

lOtiW 

Time per operation (fetch, multiply and divide) 

'i'ins ! operation 

ROM 

A^Kb 

RAM 

SKb 

Cycle/instruction basic instructions 

1 cycle 

Cycles needed for DSP multiplication 

1 cycle 

Cycles needed for integer multiplication 
Interrupt controller embedded 

1 to 4 cycles 

Other Features 

Timers, A/D converters, serial 
communication 


Table 5.1. Main Characteristics of the SH7410 (SH-DSP, SH2 Family). 
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In the light of the realities of the new CPUs, evaluating an algorithm purely by its floating¬ 


point eomplexity is elearly not enough anymore. Using a single measure is no longer possible, sinee 
the performanee of an algorithm will depend heavily on the implementation and on the arehiteeture 
of the hardware. 

The goal of the eomputational eost evaluation in this researeh is more ambitious than just 
eomparing different methods. It is intended to be used as a guideline to evaluate if a partieular CPU 
ean be used to implement the algorithm in real time. As sueh, the evaluation of eomplexity estimate 
will not measure only floating-point operations. The following items will be eounted: 


• Floating-Point Multiplication/Division fMy j . These are some of the most important oper¬ 
ations as their eost ean be mueh higher than most other operations when the eomputer does 
not have floating-point operations implemented in hardware. 

• Floating-Point Sum/Subtraction fS f). Important for the same reasons as presented above. 

• Integer Multiplication fM,-j. Integer operations ean be very fast in eomplex CPUs. But the 
low-eost and low-power ones ean usually be slower for multiplieation, taking several eyeles. 

• Integer Division fB,-). Integer divisions are a eomplex ease. Some low-end CPUs may not 
implement a full set of division operations. Thus some sort of software emulation may be 
needed. 

• Integer Sum/Subtraction (E>i). These operations are among the fastest ones, but also among 
the most frequent. For fast CPUs, they are as fast as multiplieations and divisions. But for 
low-end CPUs, summations/subtraetions are mueh faster. 

• Matrix Element Access (Matr)- Almost every operation in a eomputer oeeurs in the CPU, and 
most involve registers themselves or register and memory. Fetehing data from memory to reg¬ 
isters takes time, whieh is very important for matrix manipulations. Even taking into aeeount 
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that matrix operations frequently use the same base address, they still involve mathematieal 
operations on the index prior to fetehing the aetual value. 

• Transcendental Functions (^ranc)- The aetual implementation of these funetions varies de¬ 
pending on the algorithm used. Some dedieated hardware are ean of performing some of these 
operations, but usually they require several CPU eyeles. While this may not be pereeptible in 
fast parallelizing CPUs, the eost ean be quite high for simpler ones. Operations not provided 
by the CPU must be emulated in software and, again, the eomputation time is very dependent 
on the algorithm, implementation and desired preeision. 

C. MULTIPLE-ERROR LMS 

A referenee must be established in order to evaluate the numerieal eomplexity of the pro¬ 
posed vibration suppression method. The most natural one is the Multiple-Error LMS due to its 
widespread aeeeptanee. 

The ultimate objeetive of any vibration eontrol algorithm is to generate inputs for the plant 

//, whieh in the Multiple-Error LMS is done by using the Equation (4.4), reprodueed below: 

e-i 

W = Yj W (5- 1) 

^=0 

where x[n] is the referenee signal, Wq^i [n] is the weight matrix and Q is the order of the filter and n, 
is the input of the plant (0 < / < A, 4-1). 

The number of operations neeessary is 

r e-i 

= N< ^ 

I ?=o 

= (G-l)§/+2N{w^,, [n]x[n-^]} 

= {Q-\)Ef+Q{\Mf + mi+2Matr) 

= (G-l)§/+GM/ + G§;+2GML^r, (5.2) 
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where N{«} means number of operations required to evaluate •. 


Therefore, the number of operations required to ealeulate the whole veetor u [n] is 


N{u[n\}=Ni{Q-l)Sf+NiQMf + NiQSi+2NiQMatr. (5.3) 


The next step is to ealeulate the weights at time n+\. To do that, one of the intermediate 


steps is to evaluate rg^i for all q’s: 


N{ro,i[u]} = n| ^ 

= PS/+(P+l)N{/io,,[p]x[u-p]} 


= PSf+{P+l){lMf + lSi+2Wl,,r) 


= PSf+{P+\)Mf + {P+\)Si + 2{P+\)Wi^tr, (5.4) 


where P is the order of the estimate of the plant and ho^i [n] is the estimate of the o*^ output of the 
plant, due to the input. 


Therefore, the weight update operations eount is 


r K-i ] 

N{ Wq^i [n+\]} = N< Wqq [n] +2p Y, [n] roj [« - ^] > 

= IM^fr + 1S/+ n| 2q Y [n] roq [« - ^] | 

r K-i ] 

= IMaf;. + 1S/+ IM/ + N< Y ^o[n] ro,i [n-q] > 

= IMflfr + 1S/+ IM/ + {No — 1) §/+A^oN{eo [u] rgq [n — cf\} 


— IM^f^ + 1§/+ IMy + {No — 1) §/ + 


(IMy + !§,•+ lMa,r + N{ro,, [n — ^]}). (5.5) 


The eomputational effort to evaluate roq[n] is given by Equation (5.4) and is P^f + 
{P+ 1)My + {P+ 1)§, + 2{P+ l)Ma(r. Sinee Equation (5.5) needs time to be shifted to u — 
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one integer sum must be added to the result 


1]} = + 1S/+IMy + (A^O — 1) §y + 

A^„My + No^ i + No'Mat r + 

No (PS/+ {P + 1) My + (P+ 1) §, + 2 (P+ 1) M^fr + §,•) 

= (iVo + 1 + 2No (P+ 1)) ML,, + (1 + {No - 1) + NoP) §/+ 

{\+No + No{P+l))Mf + {No + No{P+l)+No)^i 

= (l+iVo(l + 2(P+l)))ML„+A^o(P+l)§/+ 

(l+iVo(P + 2))My+iV„(P + 2)§,-. (5.6) 

Therefore, it is possible to evaluate the eomputational effort to evaluate the matrix W (of 
dimensions QxNi) by using Equation (5.6) 

N{W[n+l]} = iV,(2(l+^o(l + 2(P+l)))ML„ + A^,(2A^o(P+l)§/+ 

NQ (1 + (P + 2)) My + NiQNo (P + 2) (5.7) 

Therefore, the total eomputational effort required to implement the Multiple-Error EMS 
algorithm is the sum of the effort to eompute u [n] and W [n + 1] 

N{W[n+l]} + N{u[n]} = (1 + (1 + 2(P+1))) ML„+iV;GiVo (P+1) §y+ 

(1 + (P + 2)) My + NiQNo (P + 3) §, + 

iV, (e - 1) §y + A^,-<2My + NQSi + INiQMatr 

= iV,e(3 + A^o(l + 2(P+l)))ML„ + 

A^,(G(l+A^o(P+l))-!)§/+ 
iV,e(2 + iVo(P + 2))My + 

A^,<2(l+A^o(P + 2))§, (5.8) 
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Note that it is safe to assume that P >> 1 and 2 >> 1 for eomplex plants. Applying these 


assumptions, an approximate estimate of the eomputational effort is 

f^{ME-LMS} ^ A,-e(3 + No(l + 2(P+l)))M«f,+ 
Ni{Q{\+No{P+l))-\)Sf+ 

NiQ{2 + No{P + 2))Mf + 

NiQ{l+No{P + 2))§i 
= 2NiNoPQMatr + NiNoPQS f+ 

NiNoPQMf+NiNoPQSi 

= NiNoPQ{2m„r + Sf+Mf + Ei). (5.9) 

The implementation of the Multiple-Error LMS used for the experiments deseribed in See- 
tion VI.D uses a plant model of order 30 and a filter order of 10. Therefore, an estimate of the 
number of operations would be 

n{ME-LMS} ra NiNoPQ{2Wl,tr + ^f+Mf + E,i) 

ra 10,800(2M«f, + §/+M/ + §,). (5.10) 

Not desiring to enter into the details of the CPU used in the experiment, it is assumed that 
the CPU takes only one eyele to perform eaeh of the eomputed measures. Therefore, the total 
number of CPU eyeles needed by the algorithm for eaeh iteration is 

n{ME- LMS} ra NiNoQ^ (2 + 1 + 1 + 1) = 54,000. (5.11) 

Sinee the sampling rate used by the algorithms in the Ultra-Quiet Platform are is IKEIz, 
that means that the CPU would need to exeeute more than 5410^ instruetions per seeond, whieh 
gives an estimate of the eloek needed. The aetual speed needed will also depend on the eompiler 
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optimization, on the code implementation and on other activities that the processor may be executing 


in parallel (I/O, D/A conversions, A/D conversions, filtering, context switching, etc). 

The microcontroller described in the previous section can execute each instruction in 33ns, 
or about 3 10^ instructions per second. Therefore, this microcontroller cannot be used to implement 
the Multiple-Error LMS on a hexapod if E > 30 and 2 > 10. 

1. Memory Requirements 

An estimate of the algorithm memory requirements is very important to decide wether a 
particular hardware setup can be used for implementation or not. For microcontrollers, this estimate 
indicates if the controller will fit in the on-chip memory or not. For more complex CPUs, this 
quantity is important to verify if the application can fit in the cache memory (much faster than 
conventional memory). 

From previous developments it is clear that the Multiple-Error FMS cannot be implemented 
in a low-power CPU. Since a full-fledged processor is needed, one wifh a large on-chip cache could 
be selecfed and ideally all fhe code and dafa would run from fhe cache and fhe speed would be 
improved by a large faclor. 

The biggesf portions of RAM memory required by fhe Mulfiple-Error FMS mefhod are 
used fo hold fhe plan! Iransfer funclion esfimafe H and fhe weighl W. Thus, fhe number of values 
fo be sfored is 


memh = {P+l)NiNo, (5.12) 

as fhe plan! has order P and dimensions A,- x No- 
Similarly, fhe weighl malrix W demands 

memw = QNi- (5.13) 
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The total memory required for a Multiple-Error LMS eontroller is approximately the sum 


of the memory needed for the estimate of the plant and weights 

raerciMELMS ~ Ni{Q+{P+ l)A/o)- (5.14) 

The number of elements to be stored for the hexapod example is approximately 

mem?3 6x (10+(30+ 1) x6) = 1,476. (5.15) 

The total memory is about 6^ytes (4 bytes/float). This is not a very large amount and 
might possibly fit in the eaehe of the CPU (if one is present). If the eode size ean be kept to a small 
reasonable amount (a few kilobytes), then the whole eode and data eould run off the eaehe and the 
speed gains would be substantial. 

D. ADAPTIVE DISTURBANCE CANCELLER 

Cost estimation for the Adaptive Disturbanee Caneeller is mueh easier than that of the 
Multiple-Error EMS approaeh. As in the previous seetion, the eomputation starts by evaluating the 
input to the plant using 

N{M[n]} = N{a[n]eos(cOc?^)+ sin(cOc^)} 

= N{a[n]eos(c0c«)} + N{6[n]sin(c0c?^)} + 1 §/—— IMy 

= 2 (2My -|- Trane + ) -|- ISy — 2W^tr ~ IMy 

= 3My -|- 2Tranc + + l§/> (5.16) 

where a [n\ and b [n] are the weights of the eontroller. 

Please note that the time needed to feteh the variables to the memory was added as a matrix 
operation, whieh is logieal for multiple-frequeney and distributed implementations. 
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The weight update equation eost ean be eomputed as 


N{a[n+1]} + N{ft[n + 1]} 


N{a [n] +qe[n] eos((Orfn)} + 

N{ [n] + qe [n\ sin(cOrfn)} — (eommon part) 

2N{a [n] + qe [n] eos((Orfn)} — AWan — 2My 
—2Mj — IMatr + 2§/ + 2N{qe[n] eos(cOdn)} 

—2Mj — IW^atr + 2§y + 2 (SMy + IT^anc + 

4My + 2T ranc + + 2E>f. (5.17) 


By summing values from Equation (5.16) and Equation (5.17), the eost is 


N{ADC} — My + 2T^a„c + 6MQf^ + 2§y+3My + 2Tra„c + + ISy 

= 4My + 4Tranc + lOM^^r + 3§/. (5.18) 

Transeendental funetions are very eostly, and it is reasonable to eompute and store them for 
use in both parts of the algorithm. Ideally, these values may be stored in registers. However, the 
worst ease, in whieh the values are stored in eonventional memory, was assumed. Thus, 


N{ADC} = 4M/ + 2T™„,+ 14ML,, + 3§/. (5.19) 

Note that N{ADC} may be halved when the proeessor used to implement the algorithm has 
a large number of registers or temporary variables are earefully arranged. 

Now, the number of eomputations required, assuming Nc eontrollers running in parallel and 
Nf frequeneies to be eontrolled, is 

N{ADC}^^^^^ = NcN/(4My + 2T™„,+ 14ML,, + 3§/). (5.20) 

It is important to realize that the above estimate did not assume the eode to be optimized, 
and the aetual eost estimate ean be improved by saving temporary results. 
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The actual count obtained for a hexapod case controlling controlling six frequencies (three 


fundamentals and three harmonics) is: 

N{ADC}6^6 = 36(4My+2T™„,+ 14M^,, + 3§/). (5.21) 

Sine/cosine functions can be evaluated by using a series expansion or a table look-up if the 
CPU does not have them implemented in hardware. The series expansion would take around 300 
cycles (?3 20 terms) if each operation could be performed in one cycle. Since the implementation 
only needs 12 such expansions per cycle, the number of operations would be 

N{ADC}6 6 36x (4+14 + 3) + 300x 12 

4,356. (5.22) 

Thus, the total time would be approximately lAA^s when the CPU mentioned in Table 5.1 
(p. 76) is used. Note that the calculation of the transcendental function evaluation is responsible for 
most of this time (82%). A table look-up would easily reduce time to a much more reasonable value 
(less than lOOq^) at the cost of added memory. 

1. Memory Requirements 

The Adaptive Disturbance Canceller controller requires very little memory, as only six vari¬ 
ables need to be stored and some of these can be shared among different controllers. The total 
number of variables needed is 


memA?,.,]v^ = 4A,-A/ -f A/ -f 1. (5.23) 

Therefore, the previous example would need 151 floating-point values, or about 604 bytes. 
Almost any microcontroller that is fast enough to implement the Adaptive Disturbance Canceller 
will have that amount of RAM embedded in the chip. 
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E. COMPARISON AND COMMENTS 


In some ways, comparing the computational and memory requirements of both algorithms 
is not reasonable, since both start from different hypotheses. The Multiple-Error LMS approach re¬ 
quires both a reference signal and a complete model of the plant transfer function whereas the Adap¬ 
tive Disturbance Canceller requires the knowledge of the disturbance frequency, no reference signal 
and an approximate model of the transfer function only at the frequencies of interest. Although 
there is some overlap, each method has different strengths. In addition, note that the Multiple-Error 
EMS method most certainly needs a full computing system to be implemented in real time, which 
means adding memory chips, interrupt controllers, A/D and D/A converters, timers, considerable 
space and much electrical power. Eurthermore, the cost increases linearly with the plant model 
and the order of the filter, and, as a result, complex systems require much more computations than 
simpler ones. 

On the other extreme is the Adaptive Disturbance Canceller, for which, in the studied exam¬ 
ple, only 4,500 operations per iteration (726 operations per assigned frequency) and 151 variables 
are needed for the six actuators controlling six frequencies each with emulated sinicos calculations. 
The Adaptive Disturbance Canceller is also very scalable: the complexity increases linearly with 
the number of frequencies to be controlled, but remains constant with the plant order. 

It is worth mentioning here that a full Adaptive Disturbance Canceller solution could be 
implemented for a hexapod with only a few chips: a microcontroller (which would include timers, 
A/D converters, D/A converters, watchdog, interrupt manager, serial communication, etc) and ana¬ 
log anti-aliasing filters. Power consumption would be quite low: most likely below IW. Einally, the 
whole controller would be small enough to be mounted in the empty space between the actuators. 
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Quantity 

Expression 

Comp. Cost 


ADC 

NcNf (4My + ITranc + 14Ma<r + 3§/) 

ME-EMS 

NiNoPQ {IMatr + §/+ M/ + §,-) 

Mem. Req. 


ADC 

4NiNf + Nf+l 

ME-EMS 

Ni{Q+{P+mo) 


Table 5.2. Computational Cost and Memory Requirement Estimates. 


Quantity 

Adaptive Disturbanee Caneeller 

Multiple-Error EMS 

Cost/iteration 

4,356 

54,000 

Memory estimate 

604 bytes 

> 6,000 bytes 


Number of Frequencies: three, one harmonic each 


Number of Inputs: six 
Number of Outputs: six 
One operation / cycle 
Floating-points implemented as floats 
Table 5.3. Adaptive Disturbanee Caneeller x Multiple-Error EMS. 
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VI. EXPERIMENTAL RESULTS 


A. OVERVIEW 

This chapter describes the implementation of the Adaptive Disturbance Canceller on the 
two hexapods described in Chapter II (p. 9). 

Section B gives an overview of some of the challenges faced when solving the vibration iso¬ 
lation problem using the Precision Pointing Hexapod. These challenges ranged from inappropriate 
design to hardware failure. 

Section C described the results on the Precision Pointing Hexapod, which exhibits strong 
nonlinearities. Experiments show that the Adaptive Disturbance Canceller algorithm and the hexa¬ 
pod configuration are not sensitive to the vibration axis. This indicates that the hexapod and the 
Adaptive Disturbance Canceller are adequate for a generic, off-the-shelf vibration isolation mount. 

The Ultra-Quiet Platform was specifically designed for vibration isolation and thus has 
better linearity properties than the Precision Pointing Hexapod, making it a better test bed for com¬ 
parison with the Multiple-Error EMS. Section D (p. 104) presents the results of this comparison. 

B. EXPERIMENTAL CHALLENGES 

This research was to be developed entirely on the Precision Pointing Hexapod, which had 
not yet been delivered when the research started. Once delivered, several tests were conducted to 
verify its behavior and also to validate the mathematical model that had been developed prior to the 
availability of the hardware. 
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The very first experiment eonsisted in exeiting the bounee mode, in whieh all the aetua- 


tors were exeited with the same sine wave. The output of eaeh aeeelerometer was then sent to 
an oseilloseope and to a speetrum analyzer. The first and most important eonelusion was that the 
system exhibited severe nonlinearity. Harmonies up to 2AKHz were observed when the input sig¬ 
nal eontained only a 50Hz sinusoidal sampled at IKHz- The experiment was repeated with a 5Hz 
sinuwave and the results were not as severe. The first experiment was relevant from the vibration- 
isolation perspeetive and the seeond experiment was important for position eontrol (pointing). The 
next experiments tested the shaker. The supplied shaker was a 25W eeramie shaker with frequeney 
response peak at 42Hz- The useful range was between 10 and lOOHz- 

Realizing potential problems with aliasing, the first experiments were made using lower 
frequeneies (10//z) in order to have all harmonies below 500Hz- Larger travel of the aetuators was 
needed at these lower frequeneies in order to obtain an aeeeleration large enough to obtain a good 
SNR. As a result, the shaker required more power and operated very elosely to its 25W speeifieation 
at the lower and of the frequeney range (lOHz)- Exeiting the shaker at sueh power levels inereased 
its own nonlinearity effeets. 

The design of the hexapod, as delivered, had the aeeelerometer as a struetural part of the 
strut, as shown in Figure 6.1. In this design, the full foree aeting on the strut passed through the 
ease of the aeeelerometer. 

While running the first experiments (low frequeney due to aliasing problems), the algorithm 
sometimes beeame unstable and the aetuators hit their limit (40A eaeh aetuator). Turning off the 
eontroller manually was not an option beeause it was implemented in software and thus it failed 
when the software entered an unresponsive state. Disabling the power supply of the aetuators man¬ 
ually was not easy beeause the switeh is loeated behind the power amplifier. Even if the aeeess were 
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Figure 6.1. Old Accelerometer Adapter. 

easier, it still depends on a fast response by the operator. Due to the original design, the force of 
the impact was transmitted through the case of the accelerometers and two of them ended up being 
damaged. 

Since two accelerometers were already damaged, all experiments were halted until the de¬ 
sign flaw was corrected. This was done by redesigning the adapters so that the accelerometers were 
inline with the strut axis but no force was transmitted though its case. The new adapter, after being 
assembled, is shown in Figure 2.7. 

Although damaging the accelerometers was no longer a concern, preventing damage to 
other components if the controller became unstable was still a concern. The range of the accelerom¬ 
eters was ±2g or 0.5 to 4.5F, assuming nominal sensitivity. It was then assumed that if the voltage 
fell outside the range of the accelerometer, then an unsafe condition occurred and preventive action 
had to be taken. The action selected was to curtail power from both the actuators and the shaker. 
In order to do that, a circuit was designed to monitor the output of one of the accelerometers. A 
window tested the voltage and cut the power supply of the actuators and the shaker if the voltage 
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was beyond safe limits. The power eould only be restored by pressing the reset switeh. A manual 


panic button was also ineluded and had the same effeet as if the aeeelerometer’s output were beyond 
aeeeptable limits. 

Although this eireuitry prevented damage to the experimental hardware, it still did not 
solve the aliasing problems. Additional hardware filters were needed. Among the several options, 
switehed eapaeitor filters were seleeted. These filters require only two eapaeitors and one inte¬ 
grated eireuit. The eutoff frequeney eould be seleeted by an external oscillator. Switched-capacitors 
anti-aliasing filters were implemented for all accelerometers. 

Once the safety device and anti-aliasing filters were implemented, several experiments were 
conducted. In order to generate acceleration levels high enough for a good SNR, the shaker had to be 
driven very closely to its power dissipation limits. To solve this problem, either the accelerometers 
or the shaker had to be changed. Therefore, the shaker was replaced by a larger model 501T. This 
reduced both the overheating and the shaker nonlinearities. 

The next important consideration was the ball joint clearance. The manufacturer stated that 
the clearance was around 2 — 4pm for each joint. From the simulations and experimental results it 
was determined that the amplitude of the vibration was in the order of 20 — 50pm. Therefore, the 
joint clearance was found to be the next problem to be addressed. The best solution was to replace 
the ball joints by flextures. Another option was to preload the joints. 

Since preloading the joints was faster and less expensive, it was the approach first tested. 
All joints were preloaded and the experiments repeated. Several experiments for which convergence 
was not achieved previously did converge after preloading the joints. This was the final setup for 
the experiments conducted in this research. 
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C. PRECISION POINTING HEXAPOD 


The hardware setup is that detailed in Seetion II.C (p. 15). The only modifieation was the 
preloading of the joints with springs, in order to reduee the baeklash present in the joints (around 
3/rm per joint, aeeording to CSA). A sinusoidal signal of 50Hz was sent in phase and with the same 
amplitude to all aetuators to evaluate the effeet of this baeklash. The power speetrum density of the 
aeeelerometers showed lines at 50Hz intervals until around lAKHz- The solution was to preload 
the joints with springs. Unfortunately, this introdueed axial forees on the aetuators axis whieh, 
aeeording to Moltran, ereated more statie frietion. 

The eontroller box was implemented exaetly as defined in the equations given in See¬ 
tion IV.C (p. 46). The parameter q was eoarsely tuned for eaeh frequeney being eontrolled. Adaptive 
Disturbanee Caneeller eontrollers were implemented in parallel and the outputs summed in order to 
eontrol more than one frequeney. 

Several different modes of vibration were tested in order to evaluate the algorithm sensitiv¬ 
ity. The definition of the modes is shown in Figure 4.12 (pg. 63). 

1. Along Z-Axis 

In this experiment the shaker was mounted direetly under the bottom plate, at the eenter. 
The vibration axis was eoineident with the z-axis. 

The first ease studied was the single tonal, at 40Hz- Harmonies due to nonlinearities were 
present on the error signal. The results obtained before and after the eontroller are shown in Fig¬ 
ure 6.2. As ean be seen, the eontroller did suppress the eontrolled frequeney {40Hz), but the vi¬ 
bration at the harmonies (80 and l20Hz) were inereased. The floor noise remained unaltered, as 
expeefed. 


93 



Accel (dB ) Accel (dB 


10 








f(Hz) 


(a) Before 



(b) After 

Figure 6.2. Z-Axis, One Tone. 
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Figure 6.3 shows the results when 40 and 80//z frequencies were controlled. As expected, 


the components of the error signal at both frequencies were reduced to the noise floor. Note that 
the second harmonic {\20Hz) increased up to about 30dB, as it was not controlled. This effect was 
specifically due to the nonlinearities of the hexapod. 

2. Tilt-Tip 

Another interesting configuration was the tilt-tip as it was simple to understand the move¬ 
ment and the displacement of the actuators. The shaker was mounted with the vibration axis parallel 
with the z-axis, with its center located on the x-y plane. This configuration did not involve transla¬ 
tional movement. 

Figure 6.4 shows results for a single tone in which the selected frequency {40Hz) was sup¬ 
pressed. Note that the static friction is believed to be responsible for the increase on the harmonics. 
As expected, there was no significant change in the noise floor. 

Next, the system was setup to control disturbances at 39Hz and 40Hz- As expected, the 
controller needed a longer time to converge. Another important finding was that the strong static 
friction exhibited by the actuators generated several harmonics multiple of IHz, which was the 
difference between the two frequencies. The two assigned frequencies were actually suppressed, 
but the energy levels on the uncontrolled harmonics were increased (Figure 6.5). 

3. Twist and Shear 

In this experiment the shaker was mounted with its center on the x-y plane, with the vi¬ 
bration axis perpendicular to it. This mainly excited the twist mode. The results, presented in 
Figures 6.6 and 6.7, were as expected. 
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Figure 6.3. Z-Axis, Two Tones. 
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Figure 6.4. Tilt/Tip, One Tone. 
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Figure 6.5. Tilt/Tip, Two Close Tones. 
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Figure 6.6. Twist, One Tone. 
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Figure 6.7. Twist, Two Tones. 
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For the last experiment, shear mode, the shaker was mounted on the x-y plane, with the 


vibration axis along the veetor x + y, where x and y are defined aeeording to Figure 2.1. This 
induees a translational movement to the hexapod (mueh more eomplex than a translation along the 
z-axis). The results shown in Figures 6.8 and 6.9 were again as expeeted. 

4. Comments 

It is very important that the algorithm performs reasonably well in most eases, in spite of 
all nonlinearities exhibited by the Preeision Pointing Hexapod. It is also important to mention that 
the sensors used in this hexapod had a very wide range ±2g, and thus the disturbanee level had to 
be quite large in order to exeite the aeeelerometers. 

Another very interesting eonelusion is that the hexapod was quite insensitive to the vibration 
type (linear/angular) and direetion. The experiments subjeeted the hexapod to disturbanees from 
several direetions. The performanee of the eontroller in all eases showed that the Stewart Platform 
is indeed a good eonfiguration for a generie vibration-isolation mount. 
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Figure 6.8. Horizontal Shear, One Tone. 
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Figure 6.9. Horizontal Shear, Two Tones. 
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D. COMPARING THE PRECISION POINTING HEXAPOD TO 


THE MULTIPLE-ERROR LMS 

It is very important to compare the results obtained with the Adaptive Disturbance Canceller 
to those obtained with the Multiple-Error LMS algorithm, and the same hardware must be used in 
order to perform a fair comparison. The Ultra-Quiet Platform, described in Section II.B (p. 10), was 
selected because it exhibits better linearity properties than the Precision Pointing Hexapod. 

1. Results 

First, an experiment was se tup with a single tone at 95Hz in order to verify the Ultra-Quiet 
Platform characteristics and how well the Adaptive Disturbance Canceller performed. The results 
of this experiment are shown in Figure 6.10. As can be seen, the harmonics were amplified, but the 
desired frequency was effectively suppressed. Figure 6.10(b) shows that the controller did not affect 
other frequencies other than the assigned one. 

Next, a second experiment was conducted in which five disfurbances (65, 95, 125, 130 and 
\95Hz) were fo be confrolled fo verify fhe behavior of fhe confroller under a more complex sifuafion. 
LMS-based fillers do nof perform well under Ihese conditions [15]. As can be seen in Figure 6.11, 
fhe Adaplive Dislurbance Canceller was able fo suppress all fhe assigned disfurbances successfully. 
If can also be seen lhal fhe error signal for fhe unconlrolled case shows frequency conlenls olher lhan 
fhe five tones. These are believed fo be a consequence of fhe cross-coupling among fhe differenl 
frequencies due fo nonlinearilies. This effecl was nof observed in fhe linear simulations and was 
observed in fhe simulalions using fhe Precision Poinling Hexapod model. 
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Figure 6.10. Adaptive Disturbance Canceller on UQP with One Tone. 
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Figure 6.11. Adaptive Disturbance Canceller on UQP with Five Tones. 

Having verified that the performance of the Adaptive Disturbance Canceller on the Ultra- 
Quiet Platform was within the expected parameters, a series of experiments were conducted in order 
to compare the Adaptive Disturbance Canceller method with the Multiple-Error LMS. 

One important characteristic of the Multiple-Error EMS is the requirement of a reference 
signal. This signal is used both to generate the estimate of the signal r[n] and to generate the 
input to the plant through the filter W. It is expected, therefore, that a noisy reference will degrade 
significantly the performance of EMS-based methods. With this information in mind the ME-EMS 
experiments were performed with both noisy and noiseless reference signals {SNR between 21 and 
15dB). 

The first set of experiments had a disturbance with a single tone at 95Hz- As seen in 

Eigure 6.12, both the Adaptive Disturbance Canceller and the noiseless Multiple-Error EMS had 
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almost the same performanee over the entire speetrum. It was expeeted that the noisy Multiple- 


Error LMS would not perform as well, and this was eonfirmed in the experiment. The noise floor 
was signifieantly raised at higher frequeneies, signaling that the Multiple-Error EMS does require a 
elean referenee signal in order to perform well. 

The seeond experiment eonsidered the sum of two elose sinusoidals (95 and 96Hz)- Ade¬ 
quately isolating elosely-spaeed tones is very important for large struetures. The beating frequeney 
is very small and struetural nonlinearities will demodulate this low-frequeney. Sinee most struetural 
modes are at low frequeneies, this ean eompromise the goal of the vibration isolation. As seen 
in Eigure 6.13, neither of the methods was able to suppress both disturbanees eompletely, and the 
Adaptive Disturbanee Caneeller had the worst performanee. One of the explanations for the poorer 
performanee of the Adaptive Disturbanee Caneeller was the tuning of the learning eoeffieient q. The 
eoeffieient was eoarsely tuned and was not optimized for any partieular disturbanee eonfiguration. 

Another important faet that the experiment uneovered was that the floor noise of the 
Multiple-Error EMS with noisy referenee was about 20 — 35dB above the noise floor of the Adaptive 
Disturbanee Caneeller and the Multiple-Error EMS with noiseless referenee. 

The last experiment set eompared both the Adaptive Disturbanee Caneeller and the noisy 
Multiple-Error EMS when eontrolling three frequeneies. As seen before, the Adaptive Disturbanee 
Caneeller introdueed several peaks in the power speetrum of the error due to eross-eoupling intro- 
dueed by nonlinearities. Eigure 6.14 shows that the Adaptive Disturbanee Caneeller did introduee 
several terms due to eross-eoupling, some about 10 — 20dB above the noise floor. The noise floor 
was also raised by 5 — lOdB eompared to the uneontrolled error. Although both methods were able 
to suppress the disturbanees, the wide-hand noise introdueed by the Multiple-Error EMS with noisy 
referenee ranged from 20 to 40dB. 


107 



2 . 


Comments 


As seen from the plots, the Adaptive Disturbanee Caneeller performanee is very similar 
to the performanee of the Multiple-Error LMS with noiseless referenee. When noise was injeeted 
in the disturbanee-eorrelated signal, the Multiple-Error EMS performanee suffered signifieantly: 
although the aetual disturbanees were suppressed, the noise floor was raised eonsiderably. The 
Adaptive Disturbanee Caneeller was not affeeted by this behavior beeause the information about the 
disturbanee (frequeney) is assumed to be known a priori and thus the referenee ean be synthetized 
without noise interferenee. 

Another very important faet was eross-eoupling due to nonlinearities among the several 
frequeneies verified in the experiment when several frequeneies were involved. Although this ef- 
feet was first notieed using the Adaptive Disturbanee Caneeller, it also oeeurred in the noiseless 
Multiple-Error EMS. The most important eonsequenee is the inability of the studied methods to 
deal with several simultaneous frequeneies (seen in the experiment shown in Eigure 6.15). This 
means that the most effeetive vibration-isolation approaeh is to isolate the vibration at its souree, 
where fewer frequeneies are present and thus aetive vibration isolation is more effeetive. 
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Figure 6.12. Comparison with ME-LMS — Single Tone. 
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Figure 6.13. Comparison with ME-LMS — Two Tones. 
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Figure 6.14. Comparison with ME-LMS — Three Tones. 
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Figure 6.15. Nonlinearity Effect on the ME-LMS — Eive Tones. 
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VII. SUMMARY AND CONCLUSIONS 


The new generation of optieal payloads imposes strieter requirements on the vibration lev¬ 
els, ehallenging the designers. These higher standards make it very diffieult to meet the requirements 
using only passive isolation and thus aetive solutions have been pursued. The Stewart Platform is 
identified in the literature as the most promising system for loeal vibration isolation, as it ean gener¬ 
ate forees and torques in any direetion. As added bonuses, hexapods exhibit very good stiffness and 
high foree-to-weight ratio. Finally, hexapods ean also be used to provide fine poinfing eapabilify fo 
fhe payload if fhe aefuafors are eapable of long sfrokes. 

Several feehniques for vibrafion isolation have been developed sinee 1981 when Widrow 
and Burgess presenfed fhe original Filfered-X LMS, buf fhe researeh of vibrafion isolation using 
hexapods inereased appreeiably only after fhe firsl half of fhe 1990’s when faster DSP ehips beeame 
available. Unforfunafely, fhe eondifions in spaee are nol as forgiving as on fhe ground: weigh!, spaee 
and eleefrieal power are sfriefly limifed. Therefore, alfhough several vibrafion isolation feehniques 
are available, fhe eompufafional resourees required limif fheir employmenf on eurrenf spaeeerafl 
designs. 


In 1998 Berfran and Monforo presenfed a very simple algorifhm eapable of performing vi¬ 
bration isolafion on SISO sysfems. Taking advanfage of fhe fael fhaf fhe hexapods used in vibrafion 
isolafion use smarf sfrufs, eaeh sfruf ean be viewed as a SISO system, and Iherefore fhis simple 
mefhod ean be employed on hexapods. The algorifhm is inexpensive enough fo be implemented 
wifh only a small fraelion of fhe eompufafional requiremenfs of fhe de facto sfandard Multiple-Error 
LMS. A deeper review of fhe liferafure revealed a laek of a sfabilify analysis for fhis simpler mefhod 
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that takes into account a generic plant. An analysis of the influence of non-controlled frequencies 


on the performance of the algorithm was also missing. 

It was not known initially if the algorithm would work properly on a hexapod due to the 
strong coupling among the struts. Therefore, it was important to simulate the hexapod in order 
to evaluate the Adaptive Disturbance Canceller in a controlled environment. Due to the unique 
characteristics of the Precision Pointing Hexapod, an appropriate model could not be found in the 
research and therefore a suitable model was derived. The model takes into account all geometric 
nonlinearities so that it can also be used for pointing. 

The Multiple-Error LMS is used a as reference throughout the work due to its acceptance 
and good performance. The Multiple-Error EMS algorithm was re-derived and the resulting equa¬ 
tions used for the computational requirements estimation. So that the Adaptive Disturbance Can¬ 
celler could be applied in spacecraft, a rigorous analysis of stability was provided and a guideline 
for the learning coefficient q selection was also included. 

Results show the method to be very selective and not to affect the error at frequencies 
unassigned to the controller. As a result, the actuators do not need to be oversized to prevent 
saturation, as is the case with the Multiple-Error EMS and other methods that use a reference signal. 

Simulations on a SISO plant were conducted to validate mathematical results and confirm 
the algorithm high selectivity. Simulations were also conducted on the Precision Pointing Hexapod 
model to verify the influence of the coupling. Results show that the coupling did not compromise 
the performance of the controller significantly. 

A full set of experiments was conducted on both hexapods. The Precision Pointing Hexa¬ 
pod, due to its ball joints, exhibited large nonlinearities. Despite the severe nonlinearities, the 
method did cancel the assigned frequencies, although the amplitude of the unassigned harmonics 
increased. 
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Experiments on the Ultra-Quiet Platform were required to eompare the results against the 


Multiple-Error EMS algorithm. The proposed method has approximately the same performanee 
of the Multiple-Error EMS algorithm when no noise is present on the referenee signal. However, 
Multiple-Error EMS performanee was signifieantly degraded, generating wide-hand noise on the 
error signal when noise is present on the referenee signal. The Adaptive Disturbanee Caneeller does 
not suffer from this problem beeause no referenee signal is used. 

The Ultra-Quiet Platform was built for vibration isolation and elose attention was paid to 
linearity. Even so, multiple-frequeney experiments showed that both methods suffered from the 
nonlinear effeets, generating several harmonies. The effeet is more dramatie in the Multiple-Error 
EMS implementation when noise is present in the referenee signal, inereasing the noise floor of the 
error signal at higher frequeneies Both methods were able to eaneel the disturbanees even in the 
presenee of nonlinearities. 

Einally, eomputational eosts for the Multiple-Error EMS and the Adaptive Disturbanee Can¬ 
eeller methods were eomputed and eompared. Traditional eost analysis uses floating-point opera¬ 
tions eount, but reeent advanees in mieroeleetronies deereased the eost of these operations and other 
operations must also be eounted. Additionally, memory requirements for both methods were evalu¬ 
ated, providing the designer with some valuable information when seleeting the hardware to be used 
for implementation. The eomputational eost required by the Adaptive Disturbanee Caneeller ean be 
smaller than one tenth of the Multiple-Error EMS requirements and similar results are obtained for 
the memory requirements. 
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A. CONTRIBUTIONS 


One of the main eontributions of this work is the approaeh from the eomputational effi- 
eieney perspeetive. Although mueh researeh has been done on vibration isolation, no work that 
foeused on eomputational effieieney was available previously. 

Researeh showed that the method seleeted laeked a stability analysis for a generie plant, 
and a stability analysis was derived. In addition, a guideline for the adaptation rate /r seleetion was 
presented. 

In addition, this work provides an analysis of erosstalk interferenee, predieting the inliuenee 
of uneontrolled frequeneies, whieh has not been addressed by the researeh eommunity until now. 

A fully nonlinear state-spaee model suitable for vibration and position simulation of large 
and small angles was developed for the Preeision Pointing Hexapod. This model shares the large 
angles eharaeteristies found in the roboties literature, without the restrietion of small aeeelerations 
and is suited for the geometry of the Preeision Pointing Hexapod. The models of the aetuators 
aeeelerometers were ineluded in the model, allowing it to be used as a tool for simulation and 
eontrol. 

Even though simulations were used in order to better understand the method, all the theo- 
retieal results were validated experimentally. The experiments on the Preeision Pointing Hexapod 
revealed the robustness of the method to nonlinearities. These experiments also reinforeed the faet 
that the method is robust to sensor noise, eonsistently aehieving suppression level below the noise 
floor. 

A performanee eomparison with the Multiple-Error EMS algorithm was eondueted, inelud¬ 
ing eases with noisy referenee, whieh substantially degraded the behavior. It was experimentally 
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shown that the Adaptive Disturbanee Caneeller delivers performanee eomparable to the Multiple- 


Error LMS in the optimal ease of noiseless referenee signal. If the Multiple-Error EMS uses a noisy 
referenee signal, its performanee is degraded and the Adaptive Disturbanee Caneeller offers better 
results. 

Einally, the eomputational requirements analysis of both the Multiple-Error EMS and the 
Adaptive Disturbanee Caneeller methods were ineluded. This analysis revealed that the simpler 
Adaptive Disturbanee Caneeller ean deliver at least the same performanee of referenee method 
using 1/10 or less of the eomputational resourees when the frequeney information is available. 


B. SUGGESTIONS FOR FUTURE RESEARCH 

Improving on the Adaptive Disturbanee Caneeller is not as easy as it might seem, sinee 
performanee improvement must be earefully balaneed against simplieity. The main researeh ob- 
jeetive was to have the simplest vibration suppressor possible, and improvements might easily add 
eomplexity that will defeat the main goal. 

The Adaptive Disturbanee Caneeller has two main weaknesses: the need for preeise knowl¬ 
edge of the frequeney and a frequeney-dependent learning rate (p). The first issue ean or eannot 
be relevant in praetieal situations. Several subsystems provide as outputs their frequeney or operate 
at a fixed, known, frequeney; ofhers do nof. If fhe frequeney is nol readily available, fhe nexf besf 
solufion would be fo measure fhe frequeney direefly, buf Ibis is nol always possible. Again. Iwo 
allernalives ean be used when fhe frequeney informalion is nol available and eannol be measured 
direefly: 1) esfimale fhe frequeney or 2) adapl fhe algorilhm fo Iraek fhe dislurbanee. The firsl 
oplion requires signiheanl eompulafional resourees fo be implemenled wilh good preeision and in- 
Iroduees delay, whieh impairs fhe response of fhe eonlroller fo ehanging dislurbanees. Eor Ihis elass 
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of problems, using an algorithm that can track the frequency by itself is highly important. This is 


one of the major areas for improvement in this method, and although much time for this research 
was dedicated to solving this particular problem, no acceptable solution was achieved. 

The other main problem is the selection of the learning rate. The algorithm, as presented in 
this research, uses a different learning rate for each frequency. Although the actual value of q is not 
critical, it is more interesting to use a single value of q for all frequencies. In order to do that, the 
transfer function of the plant at the assigned frequency must be either known or estimated online. 
This information could be used in several ways to improve the method’s characteristics. Assuming 
that the plant’s characteristics change slowly, this estimation does not need to use much computing 
power requirements, although the memory requirements will certainly increase appreciably. 

Finally, the mathematical model could also be improved. The use of quaternions would 
make its computation faster and improve its usability. Static friction should also be included in the 
model of the accelerometers, which would provide more realistic simulations at lower frequencies. 
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